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Construction and analysis of approximate models for 
electromagnetic scattering from imperfectly conducting 

scatterers 

Resume : Ce rapport est dedie a la construction ct l'analyse des conditions d'impedances 
generalises (GIBCs) utilisees en electromagnetisme comme approximations d'ordre superieur 
a une simple condition de conducteur parfait. Nous nous interessons ici au probleme 3-D 
modelise par les equations de Maxwell en regime harmonique. L'obtcntion des GIBCs se 
base sur un developpement asymptotique par rapport a la profondeur de peau a l'interieur 
de l'obstacle fortement conducteur. Ce developpement est justifie a tout ordre par des es- 
timations d'erreurs sur les series tronquees, et nous donnons des formules explicites de ses 
termes jusqu'a l'ordre 3. Ces formules sont ensuite utilisees pour obtenir les expressions des 
GIBCs. Les problcmes aux limites associes aux GIBCs sont analyses et nous etablirons des 
estimations d'erreurs relatives a ces modeles approchees en terme de la conductivitc. 

Mots-cles : GIBC, conditions aux limites effectives, forte conductivite, diffraction electro- 
magnetique, equations de Maxwell, couches limites, developpement asymptotique, modeles 
approches 
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1 Introduction 

Generalized Impedance Boundary Conditions (GIBC) have become a rather classical no- 
tion in the mathematical modeling of wave propagation phenomena (see for instance, [T^] 
and, [15j). They are used in electromagnetism for time harmonic scattering problems from 
obstacles that are partially or totally penetrable. The general idea is to replace the use of an 
"exact model" inside (the penetrable part of) the obstacle by approximate boundary condi- 
tions (also called equivalent or effective conditions). This idea is pertinent if the boundary 
condition can be easily handled numerically, for instance when it is local. The diffraction 
problem of electromagnetic waves by perfectly conducting obstacles coated with a thin layer 
of dielectric material is well suited for the notion of impedance conditions: due to the small 
(typically with respect to the wavelength) thickness of the coating, the effect of the layer on 
the exterior medium is, as a first approximation, local (see for instance, [TS], [T2], [7], 0, 

mi 

The application we consider here is the diffraction of waves by highly conducting materials in 
electromagnetism. In such a case, it is the well-known skin effect that creates a "thin layer" 
phenomenon. The high conductivity limitates the penetration of the wave to a boundary 
layer whose depth is inversely proportional to the square root of its magnitude. Then, here 
again, the effect of the obstacle is, as a first approximation, local. 

The first effective boundary conditions for highly absorbing obstacles was proposed by Leon- 
tovich. This condition "sees" only locally the tangent plane to the frontier. Later, Rytov, 
|14j . [15] proposed an extension of the Leontovitch condition, and his analysis was already 
based on the principle of asymptotic expansions with respect to the small parameter in 
the problem: the skin depth 5. More recently, Antoine-Barucq-Vernhet [2] proposed a new 
derivation of such conditions based on the technique of pseudo-differential operator expan- 
sions. However, in all these works, the rigorous mathematical justification of the resulting 
impedance conditions was not treated. 

This paper is the continuation of the work in in which we considered the case of the 
scalar wave equation. Our objective is to extend the results of to the case of 3D Maxwell's 
equations by constructing and analyzing GIBC's of order 1, 2 and 3 (with respect to the 
skin depth, the small parameter of the problem). These conditions are of impedance type 
(or H — to — E nature): they relates the tangential traces of the electric and magnetic fields 
via a local impedance operator. 

As in [11) , the construction of the approximate conditions relies on an asymptotic expansion 
of the exact solution, based on a scaling technique and a boundary layer expansion in the 
neighborhood of the boundary of the scatterer. If the organization of this paper contents is 
similar to |llj . its technicality is much higher. Moving from the scalar wave equation to the 
Maxwell system increases considerably the complexity of the problem at two levels. 
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• The first one is linked with the algebra involved in the formal construction of the 
asymptotic expansion of the exact solution (see Section 0|) . This is essentially due 
to the vectorial nature of the unknowns and the expression of the curl operator in a 
parametric coordinates system (see Section l4TTj) . The latter is based on the formulas 
proposed in |10] with some simplifications. 

• The second one is related to the mathematical analysis on the GIBCs. This is not 
only due to the fact that we have to deal with usual functional analysis difficulties 
linked to Maxwell equations (in particular trace operators and compact embedding 
properties - see Sections [5j and E] and Appendix EJ but also because we have to face 
some new difficulties in the case of the third order condition. The tangential differential 
operators that would naturally appear in the construction of the third order condition 
have not the good "sign properties" to be able to guarantee the existence of the 
approximate solution and the convergence (at optimal order) to the true solution. 
This leads us to apply various regularization procedures to construct the modified 
third order conditions (see Section l3T2|) . 

Our objectives in this work are essentially theoretical. The numerical pertinence of obtained 
conditions have already been demonstrated in [5] where, in particular, the interest of using 
a third order condition rather than a first or a second order condition is clearly shown. 

The outline of the article is as follows. Section [2] contains a description of the physical 
and mathematical diffraction problem at study with some basic stability properties of the 
solutions and asymptotic estimates with respect to the conductivity. We state the main 
results of our paper in Section [3] the GIBCs are presented in Sections 13.11 and 13.21 while 
the corresponding error estimates are given in Section 13.31 The formal construction of 
the asymptotic expansion is given in Section H] This construction is rigorously justified in 
Section [5j by proving optimal error estimates at each order. The last section is dedicated to 
the study of the boundary value problems associated with the GIBCs as well as the proof of 
optimal error estimates between these solutions and truncated asymptotic expansions. The 
main result of our paper is obtained as a combination of the results of Section [H] and Section 
[51 Some non standard technical results related to the if (curl) space (appropriate trace 
inequalities and special compact embedding properties) that may have their own interest 
have been gathered in Appendix A. 

2 Description of the physical model 

Let Qi be an open bounded domain in K 3 with connected complement, occupied by a 
homogeneous conducting medium. We denote by T the boundary of Qj and assume that this 
boundary is a C°° manifold. We are interested in computing the electromagnetic diffracted 
wave when the conductivity of the medium, denoted by cr s , is sufficiently high (5 denotes a 
small parameter). More precisely we assume that cr s — > oo as S — » and would like to study 
the asymptotic behavior of the diffracted electromagnetic field as 6 — > in order to derive 
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efficient approximate modeis to compute the diffracted waves. 

We assume that the exterior domain is homogeneous and the time and space scaies are 
chosen such that the wave speed is 1 in this medium. The eiectromagnetic wave propagation 
is therefore governed by the following Maxwell's equations: 

eW— + a s (x)E s - curlH 5 = F, in fi, 
at 



curlE 5 = 0, in Q, 




where O C K 3 denotes the propagative medium that we shall assume to be regular and 
simply connected with connected boundary (for instance an open ball), the functions o- s (x) 
and e{x) are defined by: 

(e,a 5 ){x) = 

where Sl e — f2 \ f2j and where e r > denotes the relative electric permittivity of the 
conducting medium. The right-hand side F denotes some source term that we shall assume 
to be harmonic in time : F{x,t) = Re{/(x) exp(iu;t)}, where w > denotes a given 
frequency, and where Re(z) denotes the real part of z. Hence, the solutions are also time 
harmonic: 

~E s (x, t) = Rc {E s (x) exp (iut) } , H 5 (x, t) = Re [H 5 (x) exp (iut) } , 

where the electromagnetic field {E s , H 5 ) is solution to the harmonic Maxwell system: 

( (i) {ieuj + a s )E 5 - cur\H s = /, in Q, 

{ (ii) iujH 5 + curl E & = 0, in 0. 

We assume that the support of the source term / does not touch Sl^. The system of equations 
|l| has to be complemented with a boundary condition on the exterior boundary 9f2, for 
instance we work with the following absorbing boundary condition 

E S T -H S xn = g, on dSl, (2) 

where Et '■= n x [E x n), n is a normal vector to dSl directed to the exterior of fl and g 
denoting some possible source term. 

Remark 2.1 According to (0), the boundary dSl can be seen as a physical absorbing bound- 
ary where a standard impedance condition is applied. The problem {7J 0) can also be seen 
as an approximation in a bounded domain ( namely SI) of the scattering problem in K 3 \ Sli . 
In such a case, the boundary condition on dQ has to be understood as an (low order) ap- 
proximation of the outgoing radiation condition at infinity and g a source term linked to the 
incident field. 
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Remark 2.2 In this paper, we could have treated as well the scattering problem in R 3 \ fij. 
The reader will easily be convinced that the obtained results can be extended to this case 
without any major difficulty. The only difference would lie in the reduction to a bounded 
domain. This additional difficulty is purely technical and not essential in the context of this 
paper whose main purpose is the treatment of the "interior boundary" V . 

As mentioned above we are interested in describing the asymptotic behavior of the solution 
for large a s . As suggested by the expression of the analytic solution where Qj is the half 
space, the appropriate small length parameter can be defined as: 

5:=l/Vwo* <^=> a s = 1/(luS 2 ). 

This small parameter defines the so-called skin depth: the "width" of the penetrable region 
inside the conducting medium is proportional to S. 

For the construction of approximate models in the exterior domain Q e it is useful to 
rewrite the problem UH]) as a transmission problem between (Ef,Hf) := (E s , H s )\q. and 
{E s e ,Hl) = {E S ,H%^ as follows: 

' -curlFf = /, 

iuHl + cm\E s e = 0, 
E 5 eT -H 5 e xn = g, 
E 5 e xn = Ef x n, 

' {ie r uj + ^)Ef- curl Hf = 0, in 

< iuHf + cuilE? = 0, in fii, (4) 

Hf xn = H s e xn, on T. 

We have chosen to split the two transmission conditions (namely the continuity of the 
tangential electric and magnetic fields) in such a way that the first one appears as a boundary 
condition in ^ for the interior field while the second one appears as a boundary condition 
in §3§ for the interior field. Roughly speaking, the approximate models are then obtained 
from replacing in system (j4|) the exact boundary condition on T by an approximate one, 
whose expression is derived from seeking appropriate asymptotic expansion of the solution 
in the boundary layer inside f^. 

2.1 Existence-Uniqueness-Stability 

With if (curl, O) denoting the space of functions V € L 2 (C) 3 such that curly e L 2 (U) 3 , 
where O is an open domain of M 3 , we define 

H(cax\, O) = {V € ff(curl, O) ; V T e L 2 t (dO)} (5) 



in fi e , 
in Vl e , 
on 9f2 
on r, 



(3) 
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where Vt is the tangential trace of V (cf. Section |4~T1 for more details), L\(dO) denotes the 
space of functions V G L 2 (dO) 3 such that V ■ n — on dO, where n denotes a normal to 
dO. We recall that (curl, O) is a Hilbert space with scalar product 

(WgCcurl,©) = (Wi»(0) + (curlf/, curl V) L2(0) + (U T ,V T ) L 2 (do) . 

Theorem 2.1 For given f G L 2 (Cl) 3 and g G L 2 (dCl) there exists an unique solution 
{E 5 ,H 5 ) G ff(curl,fi) x H(cml,Cl) satisfying {Hi. Moreover, th ere exists a positive con- 
stant C independent of 5 such that 

\\E 5 \\s (cullt n) + ~ 5 W&Wl'M < C (||/|| L2(0) + \\9huan)) ■ (6) 

Proof. The proof of existence and uniqueness can be found in [13] (Theorem 4.17). The 
solution E s is constructed (using the Helmholtz decomposition) as 

E s = E s + Vp s 

where p s G Hq(CI) and Eq G H (cm\,Cl) where 

H(cax\, Cl) = {V G #(curl, Cl) ; div((icr' 5 - ew 2 )F) = 0}, 

will be equipped with the norm of -ff(curl, Cl), of which it is a closed subspace. 

Taking the divergence of <JT]-i) one easily see that the function p 5 is solution to 

-Lo 2 (eVp 5 , V<t>)mn) + ^(V/, Vflvpit) = ^4>)l^) (7) 

for all cj) G Hq(CI). Choosing (j) = —p s then successively considering the real and imaginary 
parts of the resulting equality, one easily deduces that 

I|VpV( ) + ] Wvp'Wvm < c ||/|| ia(n) , (8) 

for some positive constant C independent of 5, which proves estimate ^ for the gradient 
part of the solution. It remains to get an estimate for Eq. The proof is based on the 
compact embedding of iJo(curl, Cl) into L 2 (Cl) 3 , whose proof can be found in [13], Theorem 
4.7 (which is an adaptation of the proofs in [TB] and [5])- The function Eq satisfies the 
variational formulation 

(curl curl $) L 2 (n) - w 2 (e ^, $) L 2 (n) + $)^ ( n s ) + t ,^t) l ^ 9 u) 

f 1 ( 9 ) 

= iw|(/,$) L a(n.) ~ (3,*t)l?(oo)} + uj 2 (e\7p s ,<$>) L 2 (n) - ^(V/, $) L 2 (f2t) 
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for all $ £ #0 (curl, 17), where p is the solution to ([7]). We first prove (by contradiction) 
that there exists a positive constant C independent of 5 such that 

(II/IIl^ + NIl^o))- (10) 
If not, then there would exist a sequence of data (f s ) and (g s ) such that 

ll/ 5 ||z a (n) + \\g S \\L 2 t (dn) = 1 
and such that the corresponding solution Eq satisfies 

ll-^o IU 2 (n) ->■ 00 as 6 -> 0. 
With an added tilda denoting a division by ||£ ; olU 2 (n)) one observes that 

Pollen) = 1, (11) 

and satisfies, for all <E> £ -ffo(curl, 17), 

(curl^,curl$) L 2 (0) -w 2 (£^,$) L 2 (0) + £(Eg, $)z»(n 4 ) + MM,t> ^^(dn) 
= zw{(/ 5 ,$) L 2 ( ^) - (g 5 ,$ T )L ? (ao)} +w 2 (EVp s ,$) L 2 (!i) - ^(Vp 5 ,$) L 2 (0i) . 

where is the solution to with / replaced by J" 5 . For $ = one gets 
||curl^||| 2(0) -c 2 (e^,^) i2(0) + ^r||^||i a(no rlli ?(a n) 
= z^{(/ 5 ,^) L2(0e) - (g 5 ,E 5 0T ) LUm) }+^(eVp s ,E 5 ) LHn) ~ £(Vp s , ^) L 2 (0i) . 

(13) 

Taking first the imaginary part (fT3|) and using estimate (j8]) and (fTT) . one proves that 

^\\E 5 J LHQi) +u;\\El T \\ Ll{m) <c(||/ s || L2(n) + ||5 4 || L?(an) ) (14) 

for some positive constant C independent of 6. One deduces after taking second the real 
part and using again estimate dHJ and (fTT]) that also 

||curl£ 5 || L 2 (n) < Ci + C a (||f |U= (n) + U^IL^n)) (15) 

for some different positive constants C\ and Ci independent of 5. It is then observed that 
the sequence (Eq) is bounded in 7fo(curl, 17). 

Applying (the trace) Lemma [A. II (proved in the Appendix) to Sq |n < one deduces from (fT4"|) 
and (dU) that 

WEq X n||„_i,„ as J->0. 



'H - 2 (r) 
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Also from compactness theorems in Sobolev spaces and (|15p one can assume (up to an 
extracted subsequence) that Eq x n is convergent in 2 ($£}). Hence, (the compactness) 
Lemma lA. 51 applied to the sequence Eq\q, c , shows that one can extract a subsequence (also 
denoted Eq\qJ that converges to some E weakly in i?o(curl, fi e ) and strongly in L 2 (f2 e ) 3 . 
Moreover 

E x n = onT. (16) 

Taking $ with support inside f2 e in fT^]) and taking the limit as 5 — > 0, show that E satisfies 
(both f s and g s tend to 0) 

curlcurlE 1 - lj 2 E = in f7 e (17) 
iuiEiT + curl E x n = on dCl. (18) 

The uniqueness of the solution in _ffo(curl, f2) to (|17M16I) proves that E 1 = in £l e . Since 
estimate ()14|) also proves that ||i?ol|L 2 (Oi) - * 0, we get that 

1 = \\ E Q\\L^(n) -> ||^||L2(!2 e ) 

which contradicts E = in J7 e and proves pop. 

Now take $ = Eq in (J9j) and use the imaginary part then the real part as previously done 
to deduce 

H£olli?(cu r i,Q) + ] H^IU'Cno < C (\\f\\ LHn) + \\g\\ LUdn) ) ■ (19) 
Estimate © is a straightforward consequence of fT§]) and ([5]) . □ 

2.2 Exponential interior decay of the solution 

It is shown that the norm of the solution in a domain strictly interior to f2j goes to faster 
than any power of 5. This is a first way to express how the interior solution concentrates 
near the boundary T. Let us indicate that this result will also be a consequence of the 
asymptotic analysis performed in next sections, however the methodology is more complex. 
The proof given here is a direct one and is independent of the subsequent analysis. The 
precise result is the following: 

Theorem 2.2 For any D > small enough so that SV( :— {x G Qf, B(x, v) C fii} is a 
non-empty set, where B{x 1 v) denotes the closed ball of center x and radius v, there exist 
two positive constants Cp and c p independent of 8 such that 

l|£fl|ff(curi,Qr) + ll#fllff(curi,fi.r) < c * exp(-c p / 5) { \ | / 1 1 L2 (n) + 1 1 g\ I L 2 (an) ) . 

Proof. The proof of this result follows the same lines as the scalar case treated in [TT] . 
For the reader convenience, we hereafter give the basic ideas. 
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We introduce a cut-off function </> p £ C°°(f2) such that 

4>p(x) = in f2 e , ^p(x) = /J" in f^, 
where the constant j3 v > is chosen such that 



IV^IU < j. (20) 



We set E = exp (4> p (x)/5) E s . Straightforward calculations show that 



curls' 5 = cxp ( - p (x)/(5){curlE - -V<£ p x E}. 



Hence, 



curl curls' 5 = exp ( — ipc (x) /S) |curlcurlE — ^V^p x curlE 
-i curl (V<fe x E) + iv^ 5 x (V<fe x E; 

Therefore E satisfies, after multiplying the above equality by exp(0 p /5) 
curl curl E — - (v0 p x curlE + curl (V</> p x E) 



+ [[-elo 2 + iLoa 5 )^ + ^V(t>p x (V0 P x E) j = iu)f in ft (21) 

Taking the L 2 (tt) 3 scalar product of this equation (f2"Tj) by E and using Stokes formulas 
yields, (recall that E = E 5 and 4>v = in fl e ) 

||curlE||| 2(f2i) - |{(V0 C x curlE,E) L 2 (0i) + (V<fe x E, curlE) L 2 (0i) } 
{-e r u 2 + ^)||E||| 2(ni) + £(V<fe x (V<fe x E),E)| 2(0i) 

(22) 

= ^{(/,£; 5 ) i2(nc) - {g,E 5 ) L , {m) ) - ||curl^||| 2(nc) 

+uj 2 \\E s \\ 2 L2{n ^ - iu\\E^\\ 2 L ^ {an) 

Let us denote by L$ the right hand side of the previous equality. According to Theorem 12. 1[ 
there exists a constant C independent of S such that 

|^|<C(||/||! 2(0) + ||.g|| 2 L2(ao) ). 
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On the other hand, thanks to inequality (I2U|) . and making use of the inequality \b — a\ > 
\b\ — \a\, we have the lower bound 

(- £ ^ 2 + ^)l|E||i 3(ni) + ^(V0 P x (V<fe x E), E)£ a(no 

> fl^YMk F ,A HFII 2 > 3 HFll 2 

for (5 is sufficiently small. One therefore deduces from (|22|) 

||curlE|| 2 L2(0i) + A ||E||| a(fli) < ? ||V0 p || oo ||E|| i2(Oi) ||curlE|| L2(Oi) +C(||/||i 2(o) + || 3 ||| ?(af2) ). 

Then, using again and inequality 2|a||6| < a 2 + b 2 one gets 

||curlE|| 2 L2(0i) + J-||E||| a(no < 2C (||/|| 2 L2(a) + \\g\\% {dQ) )- 

The final estimate can now be easily deduced by noticing that E s = exp ( — (3 V /<5)E in f2^, 
and by using the Maxwell equations to get estimates on H s from those on E s . □ 

3 Statement of the main results 

We shall denote by (E^' k ,H^' k ), the approximate solutions in the exterior domain Q e , the 
presence of the integer k meaning that these fields will provide an approximation of order 
0(S k+1 ) of the exact exterior electromagnetic field (E^,H^), in a sense that will be made 
precise by the error estimates (see Theorem [33}- They are obtained by solving the standard 
Maxwell equations in the exterior domain Q e 

%u)E^ h - curl H 5 e > k = f in Q e , 

iujH^ k + curlE s e ' k = in Q e , (23) 

E e,T ~ H e k xn = g on dn, 

where n denotes the normal to <9£! directed to the exterior of O, coupled with an appropriate 
GIBC on the interior boundary T of the form 

E s e ' k xn + Lj V 5 < k {H 5 * T ) = 0, (24) 

where n denotes the normal to T directed to the exterior of Q e , H S e ' k T is the tangential 
trace of H^' k , and where T> S ' k is an adequate local approximation of the H-to-E map for the 
Maxwell equations inside f2j, namely the operator: 

V s : H-?{cm\ Tl T) — ► ff^(div r ,r) 
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defined by 

V s <p = --Ef x n|r 
to 

where (Ef((p), Hf(tpf) is the solution of the interior boundary value problem 
' (ie r u + ^)Ef(<p) - curlJT/foO = 0, in fi„ 
< icoHf((p) + cmlEf(cp) = 0, in fi i; 

, Ht,A<P)=<P, on r - 

3.1 The "natural" GIBCs for k = 0, 1, 2. 

The approach that we shall use in Section 0] for the formal derivation of the GIBCs leads to 
the following expressions of V s - k (for k = 0,1,2,3), 

r V S -° = Q, 

J V 5 ' 1 = SVi, (25) 
[ V s ' 2 = SV~i + S 2 (H-C), 

where, \fi := + i^p denotes the complex square root of i with positive real part, and C 
and TL are the curvature and mean curvature tensors of T (we refer to Section 14.11 for more 
details). Note that the condition of order simply expresses the fact that the limit exterior 
problem when S goes to corresponds to the perfectly conducting boundary condition. 



3.2 The modified third order GIBC 

The same approach extended to k = 3 would suggest to take: 

V s ' 3 = V S ' 3 (26) 

where by definition (we refer to Section 14.11 for the definition of the surface operators V r , 
div r , curl r and curl r ) 

S 3 / \ 
Vq' 3 := 5Vi + 5 2 (H - C) + — ^ [C 2 - H 2 + e r uj 2 + V r div r + curl r curl r J . (27) 

However, we did not succeed in proving that such a choice was mathematically sound due 
to the presence of the second order surface operator V r div r + curl r curl r . As a self-adjoint 
operator in L 2 (r), this operator (more precisely the associated quadratic form) has no fix 
sign. This induces difficulties in the study of the forward problem via variational techniques 
and, as a consequence, the well-posedness of the corresponding boundary value problem is 
not clear: this is a new difficulty with respect to the scalar wave equation. 
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This is why we propose hereafter another third order condition, that (formally) gives the 
same order of accuracy as the one in ([25]) but admits good mathematical properties with 
respect to stability and error estimates. The reader will easily notice that the proposed 
modifications arc not the only possible ones (see for instance Remarks 13.11 and 13. 2[) . wc 
exhibit only one particular choice. We shall hereafter present the intuitive reasons that 
led us to introduce these modifications, postponing the rational justification to the error 
analysis of Section [ 



The first desirable (and probably necessary) property is the absorption property: 

Tie [ V S ' 3 ip -(pda>Q, 



for any smooth tangential vector field ip on T. Such a property is satisfied by the exact DtN 
operator and expresses the absorbing nature of the conductive medium: 



Tie [ V s ' 3 <f£dcr = -L f \Ef(<p)\ 2 dx. 
Jr ub z J n . 



It will play an essential role in proving the uniqueness of solutions. One can observe that 
this condition is satisfied by V s ' 1 and V s ' 2 . For V s ' 3 , we see that 

Tie V s < 3 = 6^— +5 2 (H-C) + —= (C 2 -H 2 + e r uj 2 ) + —= V r div r + —= curl r curl r . 
2 2y2 2\/2 2y2 

The problem comes from the operator V r div r which is negative in the L 2 sense. However, 
we can write formally 



SV2 S 3 „ 6 5 5 3 „ 

S 



(28) 



which suggests to define the real part of V 5 ' 3 as 



(1-r V r divr)" 1 +0(<5 b ), 



c x3 X 3 

Tie V 5 ' 3 = —= +6 2 (H-C) + —= (C 2 - H 2 + e r u 2 ) + —= curl r curl r 
2y2 2y2 2y2 

6 (l-S 2 VrdiVr)- 1 



(29) 



2^2 



Remark 3.1 The approximation process i28\) is analogous to the process used in the con- 
struction of absorbing boundary conditions for the wave equations, see |3 ^ for instance, 
where the Pade approximations are preferred to Taylor approximations in order to enforce 
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the stability of the resulting approximate problem. 



In 128]) . the splitting 



6V2 



2V2 2^2 



is somewhat arbitrary and could be changed into 



5V2 , s Sy/2 8V2 
(1 — a) — h a 



2 v '2 2 
for any a G ]0, 1[. Our choice corresponds to a = 1/2. 

The second modification was guided by the existence proof for the boundary value problem 
associated to the boundary condition (f24|) . We realized that it was useful that the imaginary 
part of D S ' 3 satisfies a "Garding type" inequality namely that the principal part of this op- 
erator be positive in the L 2 sense. This property is not satisfied by the imaginary part of 2?Jy 3 : 



Tm T>n' 3 — 5 



-V2 



2 v / 2 



(c 2 - n 2 



EM — 



2V2 



V r div r 



2V2 



curb 



rlr 



This time, the problem is due to the negative operator — curl r curl r . The same manipulation 
as for the real part of V s ' 3 suggests to define the imaginary part of T>f.' 3 as 



Im V 



8,3 _ 



5 S 3 
2V2 2V2 
S 



(C 2 - H 2 



■ E r ljJ — 



6 3 

2V2 



V r div r 



2V2 



(1 + 8 Z curl r curl r ) . 



Modifications (|2"5|) and (f5U]) lead us to introduce the operator 



(30) 



V 



5.3 



S^ + s 2 (H~C) + ^ r (C 2 -H 
2 2Vi 



2 ' ErLO 2 ) 



V2 
4 



5 ( (1 - 5 2 V r div r )~ 1 + 8 2 curl r curh 



[2 

+ i — 5 ( (1 + 8 2 curl r curb ) _1 - S 2 V r dh 



(31) 



which formally satisfies V s ' 3 = V S ' 3 + 0{5 5 ). 

It turns out that even if this condition is suitable for variational study of existence and 
uniqueness of the resulting boundary value problem, it did not enable us to have a direct 
proof of optimal error estimates (although we think it can be achieved by constructing the full 
asymptotic expansion associated with the associated boundary value problem) . We realized 
that the difficulties encountered in the analysis are related to the fact that the operator V s ' 3 
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is a pseudo-differential operator of order 2, whiles the exact impedance operator which maps 
continuously £f _1 / 2 (curi r , T) into i? _1 / 2 (div r ,T) is more something between an operator 
of order —1 and an operator of order 1. This gave us the idea to force our approximate 
operator to be of order by applying a regularization process (the Yosida regularization) to 
the operators curl r curl r and V r div r 



curl r curl r ~ curlr curlr (1 + S 2 curlr curlr ) 1 
V r div r ~ V r div r (1 - S 2 V r div r 



in 0{5 2 ) 
in 0{5 2 ) 



(32) 



Such an approximation is consistent with the 0(S 5 ) accuracy provided by T> S ' 3 since curl r curl r 
and V r div r are multiplied by <5 3 . Moreover, it does not affect the good sign properties of the 
real and imaginary parts of the operator since we "divide" by positive operators. Therefore, 
we propose for the third order condition the following expression 



<5,3 



U 2 + e r uj 2 ) 



V2 



S ( (1 - S 2 V r div r ) 1 + S 2 curlr curlr (1 + S 2 curh curlr 



(33) 



/2 

i - j- S ( (1 + S 2 curl r curlr ) _1 - S 2 V r div r (1 - S 2 V r div r 



Remark 3.2 The regularization process 132(1 can also be seen as an analogous to the stabi- 
lization process used in numerical methods such as stabilized finite elements or discontinuous 
Galerkin methods, in order to ensure optimal error estimates. Here again, in I3fy) . we chosed 
arbitrarily equal to 1 the regularization constant in the term in factor of S 2 . We could have 
chosen for instance 

curl r curlr ~ curlr curlr (1 + (3 S 2 curlr curh ) _1 

where (3 is a positive constant. However, the interest of choosing (3 = 1, is that we make 
appear the same operators (1 + S 2 curl r curl r ) _1 and (1 — S 2 V r div r ) _1 ) that are already 
present in the expression of T) S ' 3 ), which is of interest from the numerical point of view. 

The operator happens to have the good consistency, coercivity and continuity properties to 
lead to optimal error estimates. More precisely: 

• One can check that 

V s ' 3 = X>o' 3 + <5 5 ft 5 ' 3 (34) 
where the operator 1Z S ' 3 , given by 

(l — S 2 V r div r ) 1 (V r div r ) 2 + (l + S 2 curlr curlr ) 1 (curlr curh ) 2 ; 



7^ 3 = ^(l+z) 



RR n° 0123456789 



16 



Haddar & Joly & Nguyen 



maps continuously (T) in H^(T) and satisfies the uniform bound 

ll^llr^^HTCD) ^ L (35) 

• One can prove (see Lemma I6.1|) that T> S ' 3 has the following fundamental properties 
(obviously satisfied by V s ' 1 and V s ' 2 ) 

V^EL 2 (T), \\V 5 ' k <p\\ r < Ci 8 |M|r, Tie (V 5 ' k <p, <p) T > C 2 8 |M|£. 

with C\ and Ci strictly positive constants. These appear to be sufficient properties to 
transform the consistency properties V s ' 3 into optimal error estimates (see the proof 
of Lemma 16. 2\i . 

3.3 Existence, uniqueness and error estimates 

The natural functional spaces for the solutions of the approximate problems vary according 
to the regularity of their traces on T. We shall distinguish the case k = for which we set 

V° H = {He ff(curl, fi e ) ; (H x n )\gn € L 2 (dfl)}, V° = {E e V° H ; {Ex n)\ r = 0}, 

from the case k — 1 , 2 or 3 for which we set 

V k H = V| = H (curl, fi e ) 

(see ((5|) for the definition of _ff (curl, f2 e ))- Then we have the following central theorem, that 
uses and combines the partial results of Sections 4-6. 

Theorem 3.1 For k = 0, 1, 2 or 3, there exists 8k such that for 8 < Sk, the boundary value 
problem (Wfy, has a unique solution (E^ k , H^ k ) G V% X V%- Moreover, if B is a 

bounded domain containing fi, then there exists a constant Ck, independent of 8, such that 

\\E 5 C - E S e k \\H(cm\& c ) < Ck 8 k+1 . 

Remark 3.3 For k = 0, 1, the above theorem holds for all 8. 

4 Formal derivation of the GIBCs 
4.1 Preliminary material 

We recall in this section some well known facts about differential geometry and differential 
operators. 
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Local coordinates. Let n be the normal field denned on T and directed to the interior 
of fli . For a sufficiently small positive constant v (see condition ([3"§|) below) we define 

VtX = {x e tti ; dist(x, dn t ) < v}. 

To any x € fij' we uniquely associate the local parametric coordinates (x T ,i>) S T x (0, i/) 
through 

x = a; r + ^n, a; € f2^. (36) 

Tangential (or surface) differential operators. In what follows we deal with various 
fields defined on T: scalar fields cp (with values in C), vector fields V (with values in C 3 ) 
and matrix (or tensor) fields A (with values in £(<C 3 )). By definition: 

• A vector field V is tangential if and only if V ■ n = (as a scalar field along T). 

• A matrix field A is tangential if and only if A n = (as a vector field along T). 

For simplicity, we assume that these fields have at least C 1 regularity, but this can be 
removed by interpreting the derivatives in the sense of distributions. 

We recall that the surface gradient operator V r is defined by: 

V r ip(x r ) = V0(x r ), vV:T^R, 

where (p is the 3-D vector field defined locally in il^ by (p(xr + v n) = ip(xr)- Note that V r ip 
is a tangential vector field. We can define in the same way the surface gradient of a vector 
field as a tangential matrix field whose columns are the surface gradients of each component 
of the vector field. 

We denote by — div r the L 2 (r)- adjoint of V r : — div r maps a tangential vector field into a 
scalar field. More generally, if A(x r ) is a tangential matrix field on T, we define the operator 
AV r for a scalar field ip(x r ) by 

(AV r )ii := A(V r w). 

In the same way, we define the operator (AV r )■ acting on a tangential vector field V{x T ) 
as: 

3 

(AV r )-F:=^(AV r Fi)i, 

where the subscript i denotes the i th component of a vector in the canonical basis of R 3 . 

We then define the surface curl of a tangential vector filed V(x r ) and the surface vector curl 
of a scalar function ip(x r ) as 

curl r V := div r (V x n) and curl r ip := (V r p) x n. 
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Of course, the various operators V r , div r , 7\LV r , curl r and curl r applies in principle to 
functions denned on T. However, they can obviously be understood as (partial) differential 
operators acting on fields defined in the (3-D) domain T x (0, V). For instance, if cf> G 
C 1 (r x (0, v)), we define V r 4> G C°(T x (0, v)) 3 as follows (with obvious notation): 



W G (0,P), 



(•,")] 



ft;") 



We apply similar rules to div r , 1Z\7 T , curl r and curl r . The extension of these definition in 
the sense of distributions is also elementary (as soon as T is C°° ) . 



Geometrical tools. In what follows, and for the sake of the notation conciseness, we 
shall most of time not explicitly indicate the dependence on x r of the functions, except 
when we feel it necessary. We shall be more precise in mentioning the possible dependence 
with respect to the normal coordinate v. 

A particularly fundamental tensor field is the curvature tensor C, defined by C :— V r n. We 
recall that C is symmetric and C n = 0. We denote c%, C2 the other two eigenvalues of C 
(namely the principal curvatures) associated with tangential eigenvectors t\, t 2 (t\ ■ n = 
T2 ■ n = 0). We also introduce 

g := C!C 2 and h := -(ci + c 2 ) (37) 

which are respectively the Gaussian and mean curvatures of T, and also introduce the 
associated matrix fields: 

H = hl r and G=gl r , (38) 
where I r (x r ) denotes the projection operator on the tangent plane to T at x r . 

Let us introduce (this is the Jacobian of the transformation (x r , v) — > x - see Q36jl ) 

J(y) ( = J(v,x r )) :=det(/ + i/C) = 1 + 2vh + v 2 g, 
and we choose v sufficiently small in such a way that 

W<P, Vi r er, J(v,Xr)) = 1 + 2vh(x r ) + v 2 g(x T ) > 0. (39) 
Thus, for each v < D, there exists a tangential matrix field x r — > lZ u (x r ) such that 

(I + vC(x v )) 7lv(x r ) = I r (x r ). 

More precisely, there exists a tangential matrix field on T, A4(x r ) , such that: 

I r + vM := J(v) Tl v , Vi r er, V v < v. 
One easily sees (using for instance the eigenbasis (r%, T2, n) of C) that 

M = 2H - C and M C = Q. 
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The curl operator in local coordinates. The basic step of our forthcoming calculations 
will be to rewrite the Maxwell equations in the domain l^,by using the local coordinates. 
For this, we need the expression of the curl operator in the variables (x r ,v). It is shown in 
[TO] that the curl of a 3-D vector field V : fl" — * R 3 is given in parametric coordinates by : 



curly 



(ft„V r ) • (V x n) 



n„V r (V-n) xn-(Jl v CV)xn-d v {V xn), 



where V and V (defined on T x (0, v)) are related by 

V(x r , v) = V(x r + v n). 
This formula can be written in a more convenient form (for us), after multiplication by J{v): 

J{v) curl V = ((I + isM)V r )-(7xn)|n+ (J + vM)V r (V ■ n) 

- (C + vQ)v\ xn- J(v) d v [V x n), 
or, in an equivalent form, 

J{v) curly = (C r + v Cf f ) V - J{v) d u (V x n) (40) 
where we have introduced the notation 

' C r V = (curl r V)n + curi r (V ■ n) - CV x n 

C^V = (curl" V)n + curlf (V ■ ri) - QV x n (41) 

L curlf u := (MV r u) x n and curl" V = {MV r ) ■ (V X n). 

This expression is convenient for the asymptotic matching procedure, described hereafter, 
because we made explicit the (polynomial) dependence of the operators with respect to v. 

Functional spaces on T and trace spaces. 

We assume that the definition of H S (T) for any real s is well known. We shall denote by 

(■,-)r and (-,-) r , 

respectively the inner product in L 2 (T) 3 and the duality bracket V{Tf - V(Tf . 

Next, we introduce some notation for spaces of tangent vector fields along T. For any s < 0, 
we set: 

H°(T) = {V e H S {T) 3 I V ■ n = on T} (H?(T) = L 2 t (T) 
Ht S {T) = {V G H-%r) 3 I (V, vn) r = 0, V tp G H S (T)} (= (H^T))' 
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as well as 



fl" s (div r ,T) = {V S HI{Tf I div r V € H S {T)} 



£P(curl r ,T) = {V e H°(T) 3 I curl r V € fF(r)} 

equipped with their natural graph norms (we notice that H°(div r , T) and i?°(curl r , T) are 
often denoted by respectively iJ(div r , T) and -ff (curl r , T)). Finally, we recall the well known 
trace theorems: 

Theorem 4.1 The two trace mappings 

u e C°°(Q e ) 3 i ► u x n| r 



u e c°°(n e ) 3 



My = u — (u ■ n)n ( = n x (w x n)) 



can 6e extended as continuous and surjective linear applications from H (curl, Q e ) onto 

(div r , r) and H~ 2 (curl r , T) respectively. Moreover, ff" (div r ,T) is the dual of H~ 2 (curl 
and one /ias i/ie Green's formula: 

(curlw • u — u • curl v) dx = (u x n, vt)t = — (v x n, ur) r 
V (it,w) € if (curl, ft e ) 2 . 
4.2 The asymptotic ansatz. 

To formulate our ansatz, it is useful to introduce a cutoff function \ £ C°°(f^) such that 
X = 1 in 0^ and x = in f2j \ £l 2l/ for a sufficiently small v > 0. For this ansatz we are not 
interested in the part of the solution inside the support of (1 — x)i since we already know 
that the norm of the solution in this part exponentially decay to as S goes to (this is 
Theorem l2.2j) . For the remaining part of the solution, we postulate the following expansions: 



for x € f2 e , 
for x € r2 e , 



E s e (x) = E°(x) + S El(x) + S 2 E 2 e (x) + ■ ■ 
H*(x) = H° e (x) + S Hl(x) + S 2 H 2 (x) + ■ 
where E*, , I = 0, 1, • • • are functions defined on Q e and 

X {x)Ef{x) = Ef{x r , v/S) + S Ej(x r ,v/S) + S 2 E 2 (x r ,v/5) + ■■■ 
X (x)Hf(x) = Hf(x r ,v/S) + 6 H}{x T , v/S) + S 2 H 2 (x r , v/S) + ■■■ 
where x, x r and v are as in (|3l)|) and where Ef(x r ,n), H((x r ,rj) : T x 



(42) 



for x E fi? 
for x e nf 
K+ C satisfy 



(43) 



For a.e. x T £ T , 
For a.e. x r e T, 



+ OC 



|^(x r ,n)| 2 dn<+^ = 0, 



(44) 



\H^(x r ,r])\ 2 dr) < +00. 
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Remark 4.1 The condition {44ty imply that Ef and Hf are exponentially decreasing 
with respect to r\. 

For a.e. x r € T, lim Ef(x r ,rj) — 0, (exponentially fast) 

rj — >oo 

For a.e. x r S T, lim H((x r ,rj) = 0, (exponentially fast) 

q — >oo 

which is coherent with the existence of a boundary layer suggested by Theorem \2. 6 A 

Remark 4.2 Expansion makes sense since the local coordinates (x r ,v) can be used 
inside the support of %. 

In the next step, we shall identify the set of equations satisfied by (I?f,ifJ) and (Ef,Hf), 
£ > by writing, formally, that we want to solve the transmission problem l[5 |) -([4 ]) . 

In the sequel, it is useful to introduce the notation 

Ef{x T ,ri) := Ef(x r ,r]) + 5 El(x r ,n) + S 2 Ef(x r ,r]) + ■ ■ ■ (x r ,rj) G T x R+, 
Hf(x r ,r)) := H?(x r , 7]) + S H}{x r ,ri) + S 2 H?(x r , ry) + • • • (i r ,i|)erxM + . 

so that ansatz (|43|) has to be understood as 

x(x)Ef(x) = Ef(x r , v/5) + 0(5°°) for x e fif , 
X (x)Hi(x) = H?(x r , v/8) + 0(6°°) for x E fif . 



(45) 



(46) 



4.3 The equations for the exterior fields. 

This is the easy part of the job. The equations are directly derived from Q and we obtain 
that {El, Hi) satisfy 

iojEl - cmlHl = f k , in n e , 

iujHl +curlEl = 0, in Q e , (47) 

| T — Til x 71 — 9k, on dfl 

where we have set fo — f, go = g and fk = 0, gk — 0, for k > 1, (|4"7| being complemented 
with the interface condition 

El\ r {x r ) x n = El(x r , 0) x n, for x r e I\ (48) 

which completely defines (El, Hi) if E^(x T) 0) x n is known. 
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4.4 The equations for the interior fields. 

As indicated above, we need to compute the interior fields E\ . The principle consists into 
expressing this field in terms of the tangential boundary values of (iff), i < k by solving the 
interior equations. More precisely, we now substitute the expansion (|451 146p into the system 
(0]) and assume that the quantity: 

Hi x n = gk H e x n 

fe=0 

is known on T. We need of course to rewrite the equations of (J4j) in the local "scaled" 
coordinates 

(x r ,T) = v/S). 

Using formula (140|) with v = 5f] we obtain: 

J(S V ) (ie r uj + -L) Ef - (Cr + ^ Cf )#/ + d v H* xn = 0(5™), m T x [0, f ), 

iJ{5n)iuH* + (C r + ^ Cf - x n = 0(<5°°), in r x [0, f ), 

(49) 

These equations are complemented by the boundary condition 

Hf(x r ,0) xn + 0{5°°) =H 5 e xn (x v ), x r e T. (50) 

The substitution of |4*5l into |4*^1 [50)) leads to a sequence of problems that enable us 
to inductively determine the fields (E^,H^). The computations are relatively delicate but 
straightforward. The most difficult task is to explain the recurrence properly, which is 
the aim of this section. In Section 14.51 we shall compute explicitly the first terms of the 
expansions. 

It turns out to be very useful to make a change of unknown concerning the electric field. 
This is motivated by the observation that 

E° = 0. (51) 

This fact can be explained along the following lines: indeed from ([IS"]) and (|4T)f one deduces 
(at least formally) that: 

/• r+oo 

\\Ef\\h {ni) ~ * J J \E°(x r , V )\ 2 dr)da. 
Therefore, the a priori estimate (|10p . which says that \\Ef \\%2/q.\ — 0(5 2 ), implies E® = 0. 
The expansion for the electric field therefore starts with 6 E\ while for the magnetic field 
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Hf =^ 0. In some sense there is a natural shift of one power of S between the expansions of 
the electric and magnetic fields. This is why we introduce the "normalized" electric held: 



E< 



f 



E 



(52) 



S 11 

and we seek an expansion of the form 

if(ar r ,77) :=E°(a; r ,77) + 5E}(x r ,r)) + S 2 E 2 (x r ,v) + ■ • • (a*,*?) e T x R+ (53) 
with the correspondence 

£f +1 =E*\ k>l. 

We then rewrite (|4"5|) as a system of equations for (Ef,Hf) (we have multiplied the first 
equation by S) 

J(5rj) (ie r 5u + -) E| - (6 C r + 5 2 i] C?)Hf + J{8rj) d v H? xn = 0, in fx [0, f ), 



U{5rj) cuHf + {SC r +S 2 n Cf ) E? - J(<fy) d^E* x n = 



in T x [0, §), 



that we can rewrite by separating the li 8- independent" part, that we keep in the left hand 
side, from the remaining terms, that we put in the right hand side, as follows: 



d v H? x n 



1 



£ ^(fif^H?) in r 

2 



(54) 



- d n Et xn + iu; H? = ]T 51 A e (Ef , in T 

The linear operators {A^j , I = 1, 2, 3 } are given by: 



Aff (u, v) = C r v - 2hn (d n v x n + — 



A H (u, v) = —ie r uj u + rj C* f v — grj 2 (d v v X n H it) 

A H (u,v) = —2nh ie r ui u, 
Ah (u,v) = —2r/ 2 g ie r u> u, 



and the linear operators {A^\ £ = 1,2 } are given by: 

A^) (u, v) = — C r u + 2hr] (d^u x n — iuj v), 

Ag (it, v) = —r] Cp' f u + grj 2 (d v u x n — iuj v) . 
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Substituting (T4"5"|) and (|53p into (|5"4"|) then equating the same powers of S leads to the following 
systems: 



cLHf x n 



1 



- d n t\ xn + zu;H* = A e &t e , H*-'), in T 



(55) 



for k = 0, 1, 2, • • • , with the convention E| = Hf = for ^ < 0. 

Of course, these equations have to be complemented with the conditions (see (l50|) and (|44|)) 
fl?(ar r> 0) x n = ff e fc (a; r ,0), 



-Hoc. 



|fff( :Er ,7 ? )| 2 d7 7 <+C 



and 



|Ef (^ r ,77)| 2 d77 < +c 



(56) 



v x T g r. 

The reader can already notice how the roles of the variables r\ and x T have been separated. 
The variable x r appears as parameter for determining (E*, iff) from the previous (Ef , ijf )'s 
since, for each x r , one simply has to solve an ordinary differential system in the variable r\. 
The solutions to this inductive system of equations can be expressed in a general way using 
the result of the following technical lemma. For that purpose it is useful to introduce 

k 

P fc (r,R+;C 3 ) := {u(x r ,r 1 )=^2a j {x r ) rf , a, g C°°(r ; C 3 ) }. 

3=1 



Lemma 4.1 Let (f,g) e P*(r, R+; C 3 ) 2 and ip e C°°(r;M 3 ) ; Then the problem, 
Find (u,v) £ C OD (T;C 00 (R + )) 2 such that, 



1 



d„v x n + - u = e^" f(r), •), in Tx M+, 
to 

d v u x n + ilu v = e _v/ "' g(r), ■), in fxl + . 



with the conditions: 

V x T E r, u(x r , 0)xn = y(x r ), 



v x T s r, 



|u(a;r,?7)| 2 d?7 < +oo, / |u(^ r! ??)| 2 ^ < 
i) Jo 



(57) 



(58) 
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has a unique solution, which is of the form 

u(x r , rf) — e~^ v p(x r , rf) and v(x r , rf) = e~^ v q(x r , rf) (59) 

with (p, q) € Pfc+i(r, R + ; C 3 ) 2 and with the square root definition \Ji :— ^ (1 + %). 

Proof. This is a simple exercise on ordinary differential equations. For the uniqueness, 
assuming ip = and / = g = 0, we can eliminate v between the two equations of (I57|) and 
obtain, with ut = n x (u X n) = u — (u ■ n) n: 

d m UT — i ut = 0. 

Since we also have Ut{x t ,§) = 0, the only solution satisfying the second condition of (jS"5|) 
is ut — 0. One deduces from the second equation of (|57p that v = which in turn implies 
u = using the first equation of (|57| . 

We proceed in the same way to prove the existence of solutions of the form ([59"|) . First, the 
projection on n of the two equations directly gives 

(u ■ n) = u er^ ir i (/ • n), 

(v-n) =uje-^~ ir i (/•«). 
For the tangential components, one easily gets 



d m u T ~ i u T = St e 
d„„v T -iv T = g T e~^ irl , 



where 



' St := (d n g-V~ir)9) xn + zw/ T , e P fe (r, R+; C 3 ), 

g T :=nx (d v S - Vi ry/) - -g T , G P fc (r, R+; C 3 ). 
Thus, the key point is the resolution of the scalar differential equation: 

9^-# = s(^) e"^", inK+, (60) 

with V'(O) = -00- The unknown tp is sought in L 2 (R + ) and s(r]) is given in P m , the space of 
polynomials in 77 of degree less than m. We prove that the solution of this equation has the 
form: 

ip(r}) = e-^" 1 p{rj) with p e P m+ i. 



Indeed, after subtracting 



^0 e"^", 
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we can assume that ipo = 0. Introducing the space: 

Pn+1 ={pe P m I p(0) = 0} , 
Observing that [d m - i) {e"^ x*} = e"^ {£ x^ 2 [2V~i x + £ - l] }, we deduce that 

(d m - 1) G ^(e-^ 1 ' P ro+ll e-^" P m ). 

From the uniqueness result for L 2 solutions of (|60[) with Dirichlet condition at rj = 0, we 
deduce that the operator (9 W — i) is injective in the space e~^ v P m +i- Since e _v ^'' P OT 
and e~^ v P m +i have the same dimension, this operator is an isomorphism from e~^ v P m 
into e~^ v P m +i, which concludes the proof. □ 

As an application of this lemma we obtain the following result. 

Theorem 4.2 The fields P e fe x n G C°°(r;C 3 ) being given, there exists a unique sequence 

{(%j#)ec~(r ; c 3 ) 2 , fc = o,i,2,---} 

satisfying the sequence of problems Moreover, 

(e^"Ef, e^'Pf) G P fc (I\ M+; C 3 ) 2 . (61) 

Proof. We prove this theorem using an induction on fc. Suppose that the existence and 
uniqueness of (E^, Pf ) and the property (|61|) have been guaranteed up to k — 1 and 

(e^Ef, e V '' , i/f)eP m (r ) l+;l 3 ) 2 , m = 0, • • • , k - 1. (62) 
We also include in the recurrence the assumption (by convention P_i(r, M + ; ]R 3 ) = 0): 

' e Jlr > ( - d v %f x n + iuH™) G P m _i (r, R+; C 3 ) 2 m = 0, • • • , jfe - 1, 

j (63) 
e^" (ft,flf x n + - E>) G P m _i(r,M + ;C 3 ) 2 , m = 0,--- ,k- 1. 

According to Lemma |4~T1 to prove (|62[) and for m = k, it sufhees to show that the two 
right hand sides of |55|) are of the form: 

e -V?1 »p with pePn(r,l + ;C 3 ) 2 . 

To verify this, let us consider (u, v) satisfying: 

( e^" u, e^ 1 ' «) G P m (r,IR+;C 3 ) 2 , 

(64) 

( (- d v u x n + iwu), e^ 1 ' (S^i; x n + £ u) ) G P m _i (r, M+; C 3 ) 2 . 
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The fact that the operators A H , I = 3, 4, are polynomials of degree 2 in r), means that 

e Vir, A W g P m+2 (r,R+;C 3 ),for£ = 3, 4. (65) 

Moreover, using assumption (IM)) and the explicit form of A^j and for £ = 1,2, wc 
observe that 



( e^" Ag> (u,«), e^" 4 1} (u,v) ) e P m (r, R+; C 3 ) 2 , 
( A$ (u,v), A% ] (u,v) ) G P m+1 (r, R+; C 3 ) 2 



(66) 



Therefore, since e^" (E^.flf - *) G P fc _ £ (r, R+; R 3 ) 2 (according to (p|)). we deduce from 
([B5D that 

4 

e >/lu (E^- e ,H^- e ) S P fc _i(r,M+;C 3 ) 2 , 

;=3 

while, thanks to f()3"]) and by noticing that k — £ + 2 < k — 1 for £ > 3, we deduce from (l6"6"|) 
that 

2 



e^" ^ A£> JlM) G Pn(r,l+;C 3 ) 2 , 

2 

e^" ^ A% ] (IM,ifM) G P fe _i(r,M+;C 3 ) 2 . 



This leads to the desired property. 

To end the proof, it suffices to check the result for k — and k = 1, which will be done in 
Section [43] (see Remarks I4~3ll4~4f . □ 

4.5 Explicit computation of the interior fields for k = 1,2,3 

This section is devoted to the presentation of the technical details related to the computation 
of the asymptotic terms up to the order k — 3. In the sequel, we shall systematically use 
the following formulas, deduced from (|4"Tj) . 

(C r V)-n = curl r V T , C r V x n = curl r (V ■ n) x n - (CV x n) x n, (67) 

(CfV)-n = curl™V T) V x n = cirri? (V ■ n) x n - (QV x n) x rc. (68) 
Computation of (E° = E},H%). For = 0, ([55]) gives 

x n + - E° =0, in T x R+, 

w (69) 

- drj&l xn + iuH?, = in fx R+. 
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whose unique L 2 solution satisfying Hf T (x r ,rj) = H® T (x T ) is given by: 
t°(x r ,?/) = El(x r ,ri) = Vioj (H° X n)(x r ) e"^", 

(70) 

H*{x Tl r,)=Hl T {x T )e-^\ 
from which we deduce the useful information for the construction of the GIBCs, namely: 

E\ x n (> r , 0) = -Vi lo Hl T (x T ). (71) 

Remark 4.3 Notice that {TTj) proves in particular Theorem \4-2\ for k = 0. 
Computation of (Ej = E 2 ,U\). For fc = 1, ([55]) gives, using ([69]) 



(72) 



(73) 



^i/, 1 x n + - E- = C r Hf, in T x E+, 

UJ 

- d v El xn + iuj R\ = C r E°, in r x K+. 
We project ([72]) on n, use ([67]) and the expressions ([70]) for E° and if°, to obtain: 
E] ■ n = lo C r ■ n — lu [curl r ff° T ] (x r ) e^ 17 , in T x M+, 

F/ • n = - C r Ej • n = — ^= [curl r x n )l (a*) e - ^" in r x R+. 
W V i 

Next, we eliminate E* in ([72]) and get the following equation in El\ T 

(d 2 v - i) H\ T = nxd n [C r Hf] - - n x (C r E? x n) 

We use again ([6"T]) and ([7D|) to transform the right hand side. Using the following identity, 
that can easily be deduced from the definitions ([37]) and (155]) 

(C(V x n)) xn-CV = -2H V for all V G M 3 , 

we finally get after some easy manipulations 

(d 2 m - i) H[ T = 2ViH H° T (xr) e-^, 

whose unique L 2 solution satisfying H®(x r ,r)) — H® T (x T ) is given by: 

Hl T (x r ,ri) = (hI t { Xt ) - r, H ff e ° T (%-)) e~^\ (74) 
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Coming back to the first equation of ([72")) , we get 

Ej xn(x r ,n)=uj ( - ViHl T (x r ) + (C - H) ffJ T (i r ) + r) V~i H ff° r (a;r)) e - ^". (75) 
In particular 

E\ x n (s r ,0) = ( - V* Hl T (x T ) + {C-H) H° T (x r )) . (76) 



Remark 4.4 Notice that J7ff[ ), J7^[ ) and (f75[ ) prove in particular The 



for k = 1. 



Computation of (Ef = Ef,H?). The calculations are much harder and tedious than for 
the two previous cases. That is why we shall restrict ourselves to the main steps. Also, 
for the sake of simplicity, we shall often omit to mention the dependence of the various 
quantities we manipulate with respect to x r . 



For k = 2, |55]) gives, using ([72 



d„H? x n 



Ef 



d r ,E 2 xn + iuH 2 = r 



h- 

2 

Ei 



in T x 
in T x 



(77) 



where we have set 



H 
„2 



= C r Hi - i e r uj E? + rj (Cf - 2h C r ) J?°, 
= -C r E| -n (Cf -2hC r )E° i . 



We can go directly to the evaluation of Hf T which satisfies (apply n x d v to the first equation 
of (|77p . divide the second equation by u> and add the two results) 



(a: 



2 

7777 



i) Hi 



n x drjr^j 



E,T- 



(78) 



The next step consists in expressing the right hand side of (|78l) in terms of the previous 
(Ef , #f)'s. Using ([H71), (EH]) and the fact that ff? • n = 0, we first compute that 



nx r 2 H — nx curl r (if/ • n) — CH\ T — r\ (Q — 2hC)H® T — i e r uj (n x 
Next, we use the expressions (I70[) . ([75)1 and (I74p and the identity 
n x curl r (curl r (V x n)) = — V r (div r U) 

to obtain 



n x r H = 



1 



(V r di 



E r uj 2 )H® T — CH eT 



-Viv 



r) (3hC-g)Hl T e 



-Viri 
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J V H 



After differentiation, we get 

[-(V r div r + e r uj 2 )H° T + (ihC - Q)H° e T + JiCH\ T \ 
-77 V~i {ihC - Q)H° >T . 
In the same way, using again ((S7|) . and the fact that E° • n = 0, we calculate 



(79) 



E,T 



n x (r| x n) = -curl r (Ej • n) + (CK\) x n + rj [{G - 2/iC)E9] x n. 



Next, we notice that 

CV = —C \(V x n) x n] (and the same with 5) 

and use the expressions ([70]) . (|73| and (f75f respectively for E- x n, E, 1 ■ n and Ej x n to 
obtain 



1 



a; 



' _E.T 



curl r (curl r H° T ) + r]Vi (jp - 3hC)H° T x nj x n 
c({n - C)H° eT xn) xn-Vi (c(H^ 



x n) ) x n 



This can be written in a simplified form, using the following identities that hold for all 
Vel 3 and that are easily deduced from (f3"T|) and ((3HJ) 

' {C((H - C)V) xn)}xn= {3hC - C 2 - 2H 2 )V, 
(C(V xn))xn={C- 27i)V, 
{(mC - Q)(V xn)}xn = (3HC +Q- 6H 2 )V. 



We obtain 
1 , 



E,T 



curi r (curlr H° e T ) + (3HC - C 2 - 2n 2 )H° eT 



-y/iri 



,-y/ir] 



- V~i(C- m)H\ T e-^i + 77 V~i (3HC +Q- 6H 2 ) H° T 
Substituting ([79")) and (|5D|) in [|78p leads to the following equation 

{d 2 vn -i)H 2 T = e-^" {2V~iH Hl T + (C 2 + 2H 2 -g)H° T 

-(A r + e r u 2 ) H° e T -rjVi {m 2 - 2Q) ff° T } , 

where A r := V r div r — curl r curl r is the vectorial Laplace Beltrami operator. 



(80) 
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Since the L 2 solution to 



is given by 



ufa) = u(0) 



(5g_-i)« = (c + ^e - ^" inK+, 



-Viv 



we deduce that iff, is given by the expression 



ff? T (a r ,»?) = e-^> \Hl T r, UH\ T - JL ( C 2 - W 2 ) ff 2 T 

+ (A r + £r c 2 ) tf 2 T + £ (3W 2 - 0)2^} • 

Finally we go back to the first equation of to obtain, after lengthy calculations that we do 
not detail here, 

E? r x n = w e-*n { -V~i H 2 T + (C-H) H\ T - (C 2 - H 2 ) tf° T 



2V~i 



2yfi 

( e r u 2 + V r div r + curl r curl r ) H® T 



+ 7] (ViHHl T + i (5H 2 - 6HC + C 2 - A r - e r oj 2 ) H° e ^j 

2 ^ {m 2 -g)Hl T }. 



- V 



In particular, for r\ = 0, 



£7? T x n = w e 



5 " { - Vi ff 2 T + (C - H) H\ T -^(C 2 - H 2 ) H% T 

^-j= ( e r u) 2 + V r div r + curl r curl r ) H® T X. 

2yi ' > 



(81) 



4.6 Construction of the GIBCs 

The GIBCS of order k is obtained by considering the truncated expansions in Sl e 

E{ k :=jz^Ei and < fc := £ ^ 

1=0 e=o 

as (formal) approximations of order k + 1 of E s e and respectively (notice that k appears 
here as a subscript while it appears as an exponent in the notation of the solution of the 



RR n° 0123456789 



32 



Haddar & Joly & Nguyen 



approximate problem 
has 



). Using the "second" interface condition, namely l[15J). one 



E s Bik \ T {x r ) x n = 6 i Ef(x T ,Q) x n for x T e T. 



(82) 



1=0 



Substituting into (182|) the computed expressions of Ef(x r , 0) for I = 1, 2 and 3, respectively 
given by (jTTj) . ([76]) and (|8Tj) ) leads to an identity of the form 



and 



E s e>k x n + coV s < k [(< fe ) T ] = S k+1 <p{ on T, for fc = 0, 1, 2, 



££ 3 x n + oj V d / [« 3 )t] = S 4 <pl j0 on T, 



(83) 



(84) 



2 "I 

e,77> 



(86) 



where V S ' k , k = 0, 1,2 are given by ([25|) and X>g' 3 is given by ((2Tj) and where y>f G C°°(r) 3 , 
k = 0, 1, 2 are tangential vector fields given by 

' ¥>g =0, 

<p{ =V~iujHl T , (85) 

k ipi =V~iu H 2 e T +u(C-H) {Hl T + 6H, 
and obviously satisfy the estimates (for 5 small enough) 

||^IU|(r) < C k (s), k = 0,1,2 
where Cfc(s) is independent of 5, while ^ g C°°(r) 3 is given by 

= V~i u H\ T + lu{C-H) (ff 3 T + 5 H 2 e T ) 

+ -L(C 2 - H 2 ) (Hl T + S Hl T + S 3 iJ 3 T ) (87) 

+ — ^ ( e r oj 2 + V r div r + curl r curl r ) [H] t + S H 2 e T + 5 3 H\ t ) . 
2\/i 

The GIBC (J23J) is obtained for k = 0,1,2 by neglecting the right-hand side of For 
k = 3, the same process leads to the condition ((2"6"|) that is modified according to the process 
explained in Section [3.21 Notice that according to that construction, we have 



' ¥>3,0 



£* 3 x n + uV S ' k [« 3 )t] = <5 4 <4 on T, where <p\ = + 5 1l s > 3 [« 3 )t] , 
and using the property ([35]) of 7?. 15,3 , 

\\A { (r)<C 3 (s) (89) 

where C 3 (s) is independent of S. 
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4.7 Towards the theoretical justification of the GIBCs 

Our goal in the next two sections is to justify the GIBCs (j24|) by estimating the errors 



El - E 5 e ' k and H s e - H s e ' k , 

where (E^ k ,H^ k ) is the solution of the approximate problem f ([2"3")) . (j2~4"|) ), whose well- 
posedness will be shown in Section 16.11 (see Theorem 16.11) . It appears non trivial to work 
directly with the differences E\ - E*< k and H s e - H 5 e ' k , we shall use the truncated series 
(El k , k ) introduced in Section EP)1 as intermediate quantities. Therefore, the error anal- 
ysis is split into two steps: 

1. Estimate the differences E^ — E s e k and — H & e k ; this is done in Section O and more 



precisely in Lemma IQl and Corollary 15. 11 

2. Estimate the difference E s ek - E s e ' k and H 5 ek - H^ k ; this is done in Section |6~21 and 
more precisely in Theorem 1 1341 

Remark 4.5 Notice that step 1 of the proof is completely independent on the GIBC and 
will be valid for any integer k. Also, for k = 0, the second step is useless since 

E S,0 = E 5,0_ 



5 Error estimates for the truncated expansions 
5.1 Main results 

Let us introduce the fields E^ k (x), H^ k (x) : Q i-> C 3 such that 

k 



E 5 /(x) = i 



J2 5*Ei(x)=El 



k ■ 



t.=0 



e=o 



for x € rL 



C(x) slE i( x r, v/8) for x G «i 



H s *(x) = 



= < 



S'Hfa) = H s e>k , for x G n r: 



X (x) S £ H^x r , v/S) for x G fi< 



e=o 



where the local coordinates x r and v are defined as in Section 14.11 and the cut-off function 
X is defined as in section |4~21 These fields are good candidates to be good approximations 
of the exact fields (E s , H s ). The main result of this section is: 
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Lemma 5.1 For any integer k, there exists a constant C k independent of 8 such that 

(i) \\E S -E S / || ff(curl , n) < C k S k +K 

(ii) \\E S - E s /\\ LHni) < C k 5 k +-i, (90) 
(m) \\E S xn-E^ xn\\ H „i {r) < C k S k+1 . 

The proof of Lemma 15. 1[ postponed to Section 15. 3[ rely on a fundamental a priori 
estimates that we shall state and prove in Section 15.21 We first give a straightforward 
corollary of Lemma 15.11 

Corollary 5.1 For any integer k, there exists a constant C k independent of 8 such that: 

\ E t ~ ^f,fcllff(curl,O e ) < Ck 8 h+1 , 
\ H l ~ ffe,fclli?(curl,n e ) < C k S k+1 . 



Proof. Simply write 



El - E 5 e<k = El - E 5 e<k+1 + S k+1 E k+1 



which yields, since E 5 ek = E^ k+1 in f2 e , 

\\ E t - ^,fe||ir(curl,n e ) < ll-^f ~ E^ k+1 \\ H ^ culltn<; ) + 5 k+1 ||i?e + 1 ||.Ff(curl : Q c )- 

Using the estimate (|901i) of Lemma 15.11 we get 

\\E S e - < fc || ff(cur i,n.) < C k <S fe+ i + 5 k+1 ||£ e fc+1 || H(cui ,,n c ) < Ck 5 k+1 . 
The estimates for H s e - iff k is an immediate consequence of 

-ioj(H s e - Hl k ) + cml(E s e - Ei k ) = in Q e , 
{El - E! k ) + curl(ff e 5 -H* k )=0 in Q e . 



5.2 A fundamental a priori estimate 

The proof of lemma 15.11 relies of the following fundamental technical lemma. 
Lemma 5.2 Assume that E 5 £ ii(curl,f2) satisfies 

( curl curl E 5 - uj 2 ~E s = 0, znft e , 
{ iwEf, - curlE 5 x n = 0, on d£l, 



□ 



(91) 
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together with the following inequality 

[ (\cm\E s \ 2 - cu 2 \E s \ 2 ) dx + iu>( [ 
Jn \Jd 



< A 



^ IIE" x n\ 



\E 5 xn\ 2 ds + ^- 
+ 5 S l|E 5 |U 2(nj 



IE 5 1 2 dx 



(92) 



for some non-negative constants A and s independent of 5. Then there exists a constant C 
independent of 5 such that 



IE" 



ff(curl,n) 



< C6" 



\V 5 \\ LHni) <C5 s+ \ ||E 5 xn|| W(r) <C6 s +z, (93) 



for sufficiently small S. 



Proof. For convenience, we shall denote by C a positive constant whose value may change 
from one line to another but remains independent of S. We divide the proof in two steps. 

Step 1. We first prove by contradiction that ||E l5 || i 2( ) < C S s+1 . This is the main step of 
the proof which will use two important technical lemmas IA.1I and IA.5[ that are proven in 
the Appendix. 



Assume that the positive quantity 



X s := J-( s +!) IIE 5 ! 



L 2 (0) 



is unbounded as S — > 0. After extraction of a subsequence, still denoted E s with 8 — > 0, we 
can assume that X s -> +oo. Let E s = E s / \\E S \\ L 2 {n) (so that ||E l5 || L 2 (0) = 1). 

Our goal is to show that, up to the extraction of another subsequence, E s converges strongly 
in L 2 (n e ) and to obtain a contradiction by looking at the limit field E. 

To show this, we wish to apply to E s the compactness result of Lemma TA.5I with O = £l e . 
bmce div E s = and since E s is bounded in L 2 (f2 e ), we only need to show that: 



(i) curl E s is bounded in L 2 (fi e ), 

(ii) E s x n\gn converges in H~^(dfl), 

(iii) E s x n\r converges in H~^(T). 

We first notice that after division by HE* 5 1| x 2 (n) > the inequality (l92|) yields 



(94) 



^ | curl 



E s \ 2 — uj 2 \E 5 \ 2 



< 



dx + iu> 



\E 5 x n\ 2 ds + i 

dn o 



IE 5 1 2 dx 



A 
A* 



(95) 



E" 



\L 2 {n, 



■S~v \\E d x n\ 



H~2(r) 
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We shall now establish estimates on the two terms in the right hand side of ([9l)f in terms of 
||curlE 5 || L 2 (0) (namely inequalities |[SSJ and ifToTjl ). 

Considering the imaginary part of the left hand side of we observe that since l/A 5 is 
bounded, 

WWl^) <C&i\\V s x n\\ H _ h(T) + C 6 ||E 5 || LW 
Next, we use the trace inequality (|143p of Lemma [A. II with O = f2, to get 

HE 5 |li 2( ^ < C Si ||E 5 || J 2(ni) (||E 5 ||I 2(ni) + ||curlE 5 || J 2( ^) + C 5 ||E 5 || L2(ni) , 



which yields, after division by HE 5 ! 



L 2 (fii) ' 



|E 5 ||l 2(0i) < Ci 5 ||E 5 ||| 2(ni) +C 2 ** ||curlE 5 ||| 2(0i) . (96) 



Let K be a positive constant to be fixed later. Using Young's inequality ab < 2/3 a 3 / 2 

l£ 2 (£2,) ' 



1/3 b 3 with a = X- 1 S and & = X HE 511 5 



we get 

5||E 5 ||J 2(as) <^-l5l + ^||E 5 ||J 2(0i) . (97) 
Choosing C\K 3 = 3/2 and substituting (j9"B"|) into (j9"T|) . one deduces 

l L2(0i) <C<5l (l + ||curlE 5 ||j 2(a!) ), (98) 



IE 5 " 



which yields 

S- 1 ||E 5 || L2(0i) < C (l + ||curlE'||| 2(no ) . (99) 

Now considering the real part of the left hand side of and using the fact that HE" 5 ||i2(n) = 
1, we observe that 

||curlE 5 ||| 2(n) < C (l + <H ||E 5 x n|| fl _ i(r) + (T 1 ||E 5 || L2(0i) ) . (100) 

On the other hand, after multiplication by , the trace inequality (|143p applied to E" 5 is 
equivalent to 

<H ||E 5 xn|| H _ i(r) <c^ {r 1 ||e 5 || l2(0i) } + c {r^iE 5 !!^^)} 5 ||curiE 5 ||| 2(0i) . 

After applying the Cauchy-Schwarz inequality to the second term of the right hand side of 
the above inequality, we easily get, since 8 is bounded 

S-i ||E* x n|| ff _i < C {5^\\E s \\ L2{no + ||curlE*|| i2(0i) } , 

(101) 

< C {l + ||curlE 5 || J 2(0i) + ||curlE 5 || i2(ni) } , 
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where we used for the second inequality. Substituting (|10ip into (|100|) shows that 

||curlE s ||l 2(n) < C (l + ||curlE 5 ||| 2(0i) + ||curlE 5 || L2(0i) ) . 
This proves {MO))- We also deduce thanks to (|99|) and (fT0T|) that 

(5- 1 ||E' 5 || L 2 (0i) and <T* \\E S x n\\ H _i are bounded, (102) 

which proves in particular (|94l-(ii)) (E 5 x n converges to in (r)). This also means that 
the right hand side of remains bounded. Thus, going back to ([§!)]) shows that ||E 5 x 
n \\L 2 (dQ.) is bounded, which proves dMl-(iii)) by the compactness of the L 2 {d£l) embedding 
into 

Now we shall conclude the proof of Step 1. From ([Ml-(i)) one deduces that E 5 is a bounded 
sequence in If(curL fi), therefore, up to extracted subsequence, we can assume that E 5 
weakly converges in £T(curl, O) to some E. Considering the restriction to O e , thanks to 
(|M|) we can apply the compactness result of Lemma IA.5I and deduce that an extracted 
subsequence of E 5 , denoted again by E 5 for simplicity, strongly converges to E in L 2 (Q e ). 
On the other hand, we observe from (|102[) that E 5 strongly converges to in L 2 {VLi) 1 hence 
E = in f2j, which implies in particular 

E x n = on Y. (103) 

Passing to the weak limit in equations (|91|) one easily verify that 

curl curl E - w 2 E = 0, in fL, 

(104) 

iwE-r — curlE x n = 0, on dfl, 

The uniqueness of solutions to (|104|1 - (|103p in H (curl, f2 e ) implies that also E = in Q e . We 
therefore obtain that W converges to in L 2 (f2) which is contradiction with ||E 5 ||i2(Q) = 1. 
Consequently X s is bounded, that is to say 

||E 5 |U 2(0) <C^ +1 . (105) 

Step 2. We shall now proceed with the proof of estimates (l90|) . Considering the imaginary 
part of the left hand side of estimate (f92|) and applying Lemma [A. II (with O = f2j) yields 

\\E 5 \\h m < C (5 S+ ? ||curlE 5 ||J 2(0i) ||E 5 ||J 2(0i) + 5°+ 2 ||E 5 || L2(0i) ) . 

Using two times the Young inequality ab < l/2(a 2 + b 2 ), the first time with 

a = Si ||curlE s || J 2(ni) and b = \\E S \\ J 2(n . } , 
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and the second time with 

a= ||E 5 || L2(0i) and b = 5 s+2 , 

leads to (we also use 1 1 curl E s \\l2(qa < 1 1 curl E 5 ||l2(q) ) 

\\V S \\h( (W <C (5 2s+i + 5 s+2 (\\V s \\ LH(W +S\\cur\E s \\ L2m )). (106) 

On the other hand, considering this the real part of the left hand side of estimate (|92|) and 
using (|105|) . we get 

||curlE 5 || i2(fi) <C (s 2s+2 + 5 s+ ? \\curlV s \\l 2{(W \\E S \\ J 2(0j) + S s \\E S || £a(n<) ) , 

which gives, using Young's inequality once again, 

||curlE 5 ||i 2(0) < C (S 2s + 2 + 6 s (||E 5 || L2(0i) + 6 ||curlE 5 || i2(0) )) . (107) 
Combining (pTiH]) and (fTU?|) leads to 

\\V 5 \\h m + S 2 ||curlE 5 || 2 L2(0) < C (<P+ 4 + 8°+ 2 (||E 5 || L2(0i) + 5 ||curlE 5 || i2(n) )) , 
which yields 

||E 5 || i2( ^ )+ «S||curlE 5 || L2(0) <CS S+2 

and in particular the second inequality of (|93|) . The third inequality of ([93]) is a direct 
consequence of the first two ones and the application of Lemma I A. II in . □ 

Remark 5.1 Notice that since we simply used in the first step of the proof the fact that 
is bounded, we have proved in fact that 

Jim \\E 5 \\ L2(n) =0. 

5.3 The proof of Lemma I5TT1 

Let us introduce, for each integer fc, the error fields 

g S,k = E s _ E s,k^ H s : k = H s _ H s,k_ ( 108 ) 

The idea of the proof is to show that £ S ' k satisfies an a priori estimate of the type (|92f and 
then to use the stability lemma [5T^1 To prove such an estimate, we shall use the equations 
satisfied by {£ S ' k , H S ' k ), respectively in fij and f2 e . 

The equations in il e . It is straightforward to check that in the exterior domain fi e , the 
errors 

(£ s e ' k ,nl' k ) := (£ s ' k \n e ,H 5 ' k \n e ) 
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satisfies the homogeneous equation: 

(i) cmlH s e ' k + iuj£^ k = 0, in fl e 

(ii) curl£f' fc - iwH s e > k = 0, in Q e 



(109) 



(111) 



and 

(£ s e ' k ) T - H 5 / x n = 0, on 90. (110) 

Eliminating H 5 e > k in (fTU^I . we get 

curl (curl£,p) - uj 2 £ 5 e > k = 0, in tt e . 

curl£f> fe xn + iuj(£ 5 e > k ) T = 0, on dn. 
The equations in Qj. Now consider the restrictions to f2j and set 

(S 5 i ' k ,Hf k ) := (^ fc |o,,^ fc |oJ. 
It is also useful to introduce the fields 

k k 

El k {x r ,v) := El{x v , £), i^ fc (* r , ?? ) := £V ff^, J), 

so that using the local coordinates, we can write 

E^ k {x) = X Ef ik (x T , V ), H 5 x k (x) = X Hl k (x rtV ) in fi,-. 

Our goal is to show that (E^ k , H^' k ) satisfy the "interior equations" except that two small 
source terms appear at the right hand side, respectively due to the cut-off function \ and 
the truncation of the series at order k. We first compute that 

curLfff + iuE% k - ^E 5 x k = X (iuE^ + curLff? fe - ^4*) + V ^ X H lk, 

curl££ fc - iu;H^ k = X (cm\Ef k - iojHf k ) + Vx x E s ik . 

(H2) 

Thanks to the exponentially decaying nature of Ef k (x r ,n) and Hf k (x r ,r]) with respect to 
rj (cf. Theorem I4.2[) , the terms in factor of V% are exponentially small in 5. 

It remains to compute the terms in factor of x- These calculations are tedious, but the idea 
is simple and consists - in some sense - to do the same calculations as in Section l4~4l but in 
the reverse sense. According to (|52|) we define 

E S k-l 
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With the notation of Section 14.41 we have 

curl^fc - iuH' itk = rf k {x T , u/8) (113) 
where the function rf' k (x,?]) is given by 

2 

e=i 

Replacing Kf k and Hf k by their polynomial expansion in S, we get 

rf k = E SP M >< « - ^ A?) - ^ 3? + E 5 ' E §P A e (®l H ?) ■ 

p=0 1=1 p=0 

Using the equations (f55|) satisfied by the E^'s and Hf's, we get 

r f k = fff+EE A f<&> flf) - E E SP 4^ , #r')- 

Applying the change of index p + I — > p in the first sum, we get 

2 fc-l+^ 2 fe-1 

rf =-io;^^+E e ^4>r^rvEE^4%r^rt 

that is to say 

e=i P =k 

Paying attention to the above expression and using the form of the functions E^ and Hf (cf. 
Theorem I4.2[) . we see that 

cuilEf ifc _ iuJ Hl k = 5 k (g{ + S g 5 k>1 ) in supp X , 

where the functions <?q and g\ are of the form 

9{ q (x)=Pk, q (^,^)e- V ~ i r, ft! ,eP t (r,l+;C 3 ), 3 = 0,1. (114) 

From (|114p . we easily deduce that 

\\X9lj&{fU)<Ck, q 6l, ||X curl^Jl^o,) <C' k<q Si, q = 0,l. (115) 
In the same way, using again local coordinates, we have 

iuE s i k + cm\H s i k - ^El k = - sf k (x, v/S), 
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with 



4 k = dnHtk x n 



Replacing ¥,f k and Hf k by their polynomial expansion in 8 we get 



i , ft 
fc-1 



fc-1 



8,k 



p=0 



1=1 p=0 



S k d v Hf x n - X ^ 4/(0, A?) 



Using equations ([55)1 satisfied by the E*"s and Hf's, we get 

4 fc-1 4 fc-1 



<5,fc 



E E ^(Er^rVE E ^^(^n 



= 1 p=0 



= 1 p=0 



or equivalently 



<S fe d v Hf x n - E ^ +fe 



s f = "e 1 s p E^c^r^n 
P =fc ^=i 

£=i 



This time, we see that we can write 



iujEf 



1 1 3 

\ k + curLff^ - -El k = - J] 6" h s k q in supp X , 



9 =o 



where the expression of h k q is similar to the 3fc. ? 's (see formula (JTT4J) and implies in par- 
ticular that 

11x^11^(0.) <C k , g SK g = 0,l,2,3. (116) 
In summary, taking the difference between (|112p and ((¥]) we have shown that 

3 



q=0 



Vx x H[ k 



curie,- — ILOHa 



± 



q=0 
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where eliminating 7if k we get 



curl curlf i - u) ^ + — £ .' = f k 



(117) 



with 



ft ■■= * 



1 1 

k x(Y, Sq c ™ l d, q ) + Vx x (J2 6< g{ q ) + V X x curl (V X x E[ 



9=0 



9=0 



ILO 5 



k- 



9=0 



Taking into account the form of the functions gf. and the exponential decay of the the fields 
Ef and Hf fTheorem l4.2j) . and since the support of V X is separated from T, there exists a 
constant r > such that: 

V X x curl (V X x El k ) || L2( n e) < C x {k) e^ 5 , 
i 

V X x (22 S*,,) IU'(n.) < C 2 (fc)e^ 5 , 

9=0 

|V X x^ fc || L2( o e) <C 3 (fc) e-^. 
Combining these inequalities with estimates (|115|) and () 1 16[) . we see that: 

||/ fe 5 |U=(o e ) < C k 6 k -i . (118) 

Error estimates. We can now proceed with the final step of the proof. First, we multiply 

the equation (jllip by £^ k and integrate over f2 e . Using the Stokes formula and the boundary 
condition in (jllip . we get 



|curl£f< fc | 2 dx 



£ 5 > k \ 2 dx 



Next, we multiply the equation (|117[) by £ - ' and integrate over Oj. We get 



i,fc 



x n| 2 + (curl£ e 5 > fe x n, (£ 5 e ' k ) T ) i 



= 0. 



|curl£p| 2 dx-u 2 



o, 



£^dx-i-l |^ fc | 2 -((curl^xn).(^ fc ) 



1 



S2, 



a, 



~S.k\ 



ft-£t' k dx. 



Adding the last two equalities and using the fact that £ s,k belongs to H (curl; ft), we get 



/ |curl^ fc | 2 -^ 2 / \£ s > k \ 2 -iJ [ \£ S ' k xn\ 2 + ± [ \£ s < k \ 2 



cml£ s e > k x n - curl£*' fc x n, (£ s > k ) T ) + I f k ■ £-' k dx 



(119) 
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It remains to compute the jump 

curl£f' fe X n- cm\£f k Xn= cm\E 5 ek x n - curl£f jfc x n 

across T. Taking the trace on T of equation (|113p . we get, with p^ fe (a; r ) = rf k (x r , 0), 

curlE? fc x n = iujH- k xn + pf k on T. 
The function p S { ' k is not zero but small. In particular, according to , we have 

\\pf k \\ Hi{r) <C k 5 k . (120) 

On the other hand, taking the trace on F of the first equation of (|109|) we get 

curLEf k x n = iujH^ k x n, on V 



The continuity conditions (I56|) imply H ik x n = H^ k x n on T so that 

curl£f' fe x n - curl£f fe x n = cm\E s ek xn- curlE^ x n = pf k . (121) 
Substituting (|12ip into (|119p and using the estimates (|118p and (|120p . we get 

/ |curl£ 5 > fc | 2 - c 2 / \£ s < k \ 2 -iJ( \£ s * x n\ 2 + 1 / |^< fc | 2 ) 

< C fe ( 5 fc "i x n\\ L2((W + S k \\£ 5 > k x n\\ H _ i{r) ) . 

We can finally apply Lemma I5T21 with E s = £ S ' k and s — k — |, which provides the desired 
estimates. □ 

6 Analysis of the GIBCs 

6.1 Well-posedness of the approximate problems 

We shall prove in this section that the approximate fields [E °' , H ' ), solution of (f2"31 12~4"]) 
for k = 0, 1, 2, 3, are well defined. In fact, for k < 2 this result is an application (or 
an adaptation) of classical results about Maxwell equations with an impedance boundary 
condition of the form 

E xn + uZ H T = onT, 

where Z is a function with positive real part (see for instance |13|). To include the case 
k = 3 it is sufficient to extend these results to the cases where Z is a continuous operator 
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form L\ (r) into L\ (T) with positive definite real part. More precisely we shall assume that 
there exists two positive constants and z* such that 

f (<) \\Z<p\\ T < z* Mr, 

{ (122) 

for all 6 ^t(r)- These properties are satisfied by the operators T> S ' k , k — 1,2,3 for 
S sufficiently small, and can be seen as a special consequence of Lemma 16.11 (stated and 
proved in next section) where the dependence of the constants and z* in terms of 5 is 
also given (which is important for error analysis). The functional space adapted to this type 
of boundary conditions is the same as for constant impedances, namely H (curl, f2 e ) (see ([5]) 
for the definition of this space). 

Theorem 6.1 Let f € L 2 (f2 e ) be compactly supported in fl e and g £ L^(dfl e ). Then the 
boundary value problem 

curli/ + iloE = 0, in 51 e , 
curli? — iujH = 0, in f2 e , 
E x n + Ht = g, on <9f2, 

E x n + ujZHt = .<?, on T. 
has a unique solution (E,H) in lf(curl, Q e ) x iJ(curl, f2 e ). 

Proof. The proof uses basically the same arguments as for classical impedance conditions 
(see for instance [O]) and is detailed here for the reader convenience. The approach is 
divided in three steps: 

1 . One eliminates the electric field E and formulate a boundary value problem in H only, 
then writes the associated weak formulation in an appropriate functional framework. 

2. One shows that the Fredholm alternative can be applied for solving this problem. 

3. One shows the uniqueness of the solution which also implies the existence. 
Step 1. The problem to be solved for H £ if (curl, fi e ) is 

curl curl if — u> 2 H — —curl/ in fi e , 

curlH xn-i lu 2 ZH t = -iwg, on T, (123) 
cuvlH x n — i lu Ht = —i^g, on <9f2. 
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The equivalent weak formulation of the boundary value problem (II 23[) can be written in the 
form: 

Find if £ if (curl, f2 e ) such that, 

J (curlif • curl if' - lo 2 H-TT) dx - % uo{H T ,H' T ) dn - i u 2 {ZH T ,H' T ) T (m) 

= -iw{g,H! r ) me + [ /-curlif' dx, V if' G H (curl, fi e ). 

Next, we formulate a variational problem, equivalent to (|124[) , posed in a subspace ffo(curl, f2 e ) 
of if (curl, fi e ) having additional compactness properties, by using the Helmholtz decompo- 
sition. Let us consider the closed subspace of if 1 (0), 

S:={ve H\n) I («| r ,ti|«i) g P (r) x P (dn)} 

where f"b(r) is the space of constant functions on T (the same for Po(dft)). Then we define 
ffo(curl, f2 e ) as 

ff (curl, fi e ) :={ffG ff (curl, fi e ) / [ if • Vw dx = 0, V v G 5}, 

which forms a closed subset of if (curl, f2 e ). One also has the orthogonal decomposition (in 
if (curl, fl e ) = Ho (curl, ft e ) © V5, where VS := {Vw / v G 5}. (125) 



Remark 6.1 By Green's theorem, the reader will notice that ffo(curl, Cl e ) is nothing but the 
subspace of if (curl, fi e ) made of vector fields whose divergence vanishes and whose normal 
trace on each connected component of dfl e , namely T and dfl, has zero mean value. 

We claim that the problem Q124p is equivalent to 
Find if G ifo(curl, f2 e ) such that, 

/ (curlif • cmlH 7 - u?R -if 7 ) dx - i oj(H T ,H' T ) dn - i oj 2 {ZH T ,H' T ) T , , 

= - i u{g,H' T ) d n e + [ /-curlif' dx, V if' G H (curl, fi e ). 

Indeed, let ff G ff (curl, fi e ) be a solution of (j 1 24|) . Since, ffo(curl, f2 e ) C if (curl, fi e ), to 
prove that ff is solution of (|126j) . we only have to prove that ff G ffo(curl, f2 e ). According 
to (|125|) . we have: 

H = H + Vu, if G ff (curl, fi e ), «eS, (127) 
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and we simply have to prove that Vu = 0. If we choose H' = Vtt with v G S in (|124p . 
remarking that curlii' = and H' T = on Y and dtt, we get, using (|127[l : 



|Vu| 2 dx = =4> Vu = 0. 



Reciprocally, one proves that any solution of (|126p is solution of (|124[) using the same type 
of argument (decompose the test function instead of the solution) . 

Step 2. We now rewrite (|126|) in the form 

{Find H £ i?o(curl, f2 e ) such that, 
(128) 
a{H, H') = b(H, H') + L(H'), V if' € # (curl, fi e ), 

where we have defined: 

a(H, H') = I ( curl H ■ curl W + H ■ TP) dx - i uj{H Ti H' T ) dn - i oj 2 (ZH T , H' T ), 

b(H, H') = (1+w 2 ) / H W dx, 
Jn c 

L(H') = - i uj(g, H^) d n e + [ f ■ curLH 7 dx. 

According to Q122[) , one has, with a < 1 denoting an arbitrary positive constant 
\a(H,H)\ > %\Kea{H,H)\ + \\lma{H,H)\ 

> l(\\H\\ 2 H(curW - z*\\Ht\\ 2 t) + ^\\Ht\\ 2 t + mrWln- 

Choosing a such that az* < uj 2 z* proves that the sesquilinear form a is coercive on 
i?o(curl, Q e ). On the other hand the sesquilinear form b is continuous with respect to the 
L 2 (fl e ) norm, which compactly contains 7Jo(curl, fi e ). Therefore the Fredholm alternative 
can be applied to (|128p : existence uniqueness and stability is equivalent to the uniqueness 
of solutions. 

Step 3. We shall now prove the uniqueness of the solution of (1124p . Let H be a solution of 
(|124p with / = and g — 0. Let us take H' — H and condider the imaginary part of the 
resulting equality. We get 

lu \H\ 2 dn + lu 2 Tie (ZH T , H t ) = 0. 

From (|122p one deduces that Ht = on dft and, using the boundary condition on dfl that 
curliJ X n = on dil (this makes sense since, curl cm\H = uj 2 H, and therefore curl if 
belongs to -ff(curl, f2 e )). One concludes using standard unique continuation theorems for 
Maxwell's equations (see [13]). □ 
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6.2 Error estimates for the GIBCs 

The error estimates rely of some key properties of the boundary operator V S ' k that we shall 
summarize in the following lemma. We recall that V 5 ' k = 0, Sy/i and 6yi + 5 2 (Tt — C) for 
k — 0, 1 and 2, respectively. For k — 3, we denote by A s and B s the two operators 

A 5 := (l - 5 2 V r divr B s = (l + 8 2 curl r curl r ) _1 . 

By Lax-Milgram's Lemma these operators are well defined as continuous operators from 
L 2 (r) to respectively H(div r , T) and if(curl r , T). Setting 

a 5 := -L + 5(H - C) + ^= (C 2 - H 2 + e r u, 2 ) , 



1 6 



2 



the expression (|3"3"|) of I?" 5,3 can be written in the form 

V 5 < 3 cp = 6a S Lp+^- 6 (A s (f + S 2 cml r curl r B 5 ip 



/2 

+ itfftp + i-^- 5 (SV - <5 2 V r div r ^V) ■ 

The fundamental properties of the operators T> s,k are summarized in the following lemma. 

Lemma 6.1 Let k = 1,2 or 3. There exist a constant 5k > and iwo constants C± > 
and C2 > 0, independent of 5, such that 

(i) ||X*Vllr < Ci * |M|r, 

129 

(m) fte (r> 6 'V^)r > C 2 J ll^lll, 
/or a/Z </? € Lj(r) and <5 < <5fc- 

Proof. These properties are straightforward for k — 1 and 2. We shall concentrate on the 
case k = 3. We first observe that a 5 and (3s are bounded functions on T, and if 99 € £?(r) 
then <S 2 V r div r AV = (^V - </>) G ^t(r) and <5 2 curl r curl r B<V = {-B s <p + <p) e L 2 (T). 
Therefore P 5 ' 3 ^ S L 2 (T) and one has 

(X>*'V,V0r = *(o tf ^,V')r + ^*((^Vi^)r + « 2 (curlr curl r BVV'Jr) 

(130) 

+ l( S (/3 a v», V)r + *^ <* ((^V VOr - <S 2 (V r div r A s <p, ^) r ) 
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for all tp, tp G i|(r). For <$ sufficiently small, the functions 0:5 and f3& satisfy 



< a* < \a$\ < a* and < /3» < \(3 S \ < (3* 



(131) 



for some positive constants a*, a*, /3* and /?* independent of <5. On the other hand, from 
the identities 

(1 - <S 2 V r div r )A s tp = tp and (1 + £ 2 curl r curl r )B s tp = tp 
one respectively deduces 



and 



(A s tp., tp) r 
-(V r div r V)r 

(£?W)r 
(curl r curl r B s tp, tp)p 



= \\A^f r + S 2 \\diy r A s tpf r 
= ||div r ylVllr + ^ 2 ||Vrdiv r ylVllr 
= \\B 5 tp\\^ + S 2 \\cuxl T B 5 tp\\^ 



(132) 



(133) 



1 1 curl r B s tp\\ 2 T + (5 2 ||curl r curl r B s tp\\l 

Property (it) is obtained as an immediate consequence of (|132[) . (|133[) and (| 1 3 1 [) when applied 
to (|130[) with tp = tp. Identities (|132[) and (|133[) also respectively imply, 

\\A s tp\\ r < \\tp\\ r , <5 2 ||V r div r A s tp\\ T < |M|r, 

\\B 5 tp\\r < llvllr and 5 2 \\cml r curl r B 5 tp\\ r < \\tp\\ r . 

Property (i) is then easily obtained from p30|) with tp = T> S ' 3 tp and using these estimates, 
as well as ([TgTjl . □ 



6.3 Error estimates for the GIBCs 

We shall set for k = 0,1,2,3, 

cS,k j?5,k y-vfc j^t 

°e — 1=0 "^e; 

-ljS,k _ TjS,k _ rrf 



(134) 



Using ([83]), together with $M§ and ((35]) when fc = 3, we see that (£ s e k ,H s e k ) € V| x is 
solution of the boundary value problem: 



curl H 5 e > k + iuj E^ k = 0, 
curl £ s - k - iuj U s - k = 0, 



Hf< fc x n = 0, 



7l5,fc 



x n + ojV s > k (H s e ' k ) T = S k+1 tp{ 



in fi e , 
in fi e , 
on c?fl, 
on r, 



(135) 
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where the tangential vector fields ipt remain bounded with respect to S in all spaces i/ t s (T) 3 . 
Eliminating £*' , we see that H s e < k € V k H satisfies 

' curl(curl7if< fc )- w2 ^e fc = in fl e , 

< curltt*'* xn-iuj 2 V S ' k {H s e ' k ) T = 6 k+1 <p s k , on T, 

K curl ftp xn-iw (W*' fc ) T = 0, on 90. 

The proof of error estimates is based on some key a priori estimates that we shall give 

hereafter. We multiply by ftf' fe the equation satisfied by ftf' fe in f2 e , integrate over O e and 
use Green's formula to obtain, after having used the boundary conditions on d£l and T: 

I (|curlftf \ 2 -u 2 \H 5 e < k \ 2 )dx-i [ \(H s - k ) T \ 2 da 

Ju e Jan (136) 

- i co 2 {V s ' k (H s e > k ) T , (ftf' fc ) T ) r = S k+1 (y k , (Wp) T ) r , 

where ( , ) r here denotes a duality pairing between i/~ 1 / 2 (div, T) and i/~ 1//2 (curl, T). Con- 
sidering the imaginary part of (|136p and using l|129j) -(ii) together with trace theorems in 
H (curl, Q e ), one obtains the existence of two non negative constants G\ and Ci independent 
of 6 such that 

Cx 8 \\{H^ k ) T \\ 2 + \\(H 5 e > k ) T \\ 2 dn < C 2 5 k+l \\H 5 > k \\ H{cwX ,n B ). (137) 

More precisely we have C\ = for k = and Ci > for k ^ 0. Using (fH9"|) -fi) and (fl37)) 
one also deduces that 

\Im {V s > k (Hi' k ) T , (H s e > k ) T ) r \ < C 3 S k+1 ||Wp||H(cun,o e ), (138) 

for some constant C3 independent of 6. Now considering the real part of (|136p and using 
(|138|) as well as the trace theorem in ft(curl, fi e ), one gets the existence of two positive 
constants C4 and C5 independent of <5 such that 

ll^f' fe Hff(curl,n e ) < C4 5 k+1 ||ftf' fc || ff (curl,O c ) + C5 IIKPllI^)- (139) 

Based on these a priori estimates we are in position to prove the following result. 
Lemma 6.2 For k = 0, 1, 2 or 5, £/iere eiisi a constant C independent of 8 and 5q > swc/i 

ll^f' fc ||H(curl,n e ) + ||We fe ||ir(curl,fi e ) < Cfc 8 k+1 

for all 5 < do- 
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Proof. According to estimate (|139[) and first two equations of (|f 35|) it is sufficient to prove 
the existence of a constant C independent of 8 such that 

\\n 5 A\ LH n c) <C5 k+l . (140) 

Let us assume that (|140[) does not hold, i.e. Xs := £ fc+1 /ll^e fc |U 2 (n e ) g° es to as S — > 0, 
and consider the scaled fields 

fc a = WP/l|WPll^(Oe) and = £ 5 e > k /\\H 5 e > k \\ L 2 { n e y 

Dividing (|139p by H^f' fe |||2(n e ) implies in particular that (h 5 ) is a bounded sequence in 
i?(curl, f2 e ). Dividing (| 1 37[) by the same quantity and using the latter result shows that 

Ci 5 \\h s T \\l + \\h 5 T \\ 2 dn -> as S -► 0. (141) 

The last boundary condition in ()135j) combined with property (|129|) -(i) shows in particular 
that 

|l4||r<C 6 5||/4||r + A 5 ||^|| r 

(with the alternative, Cq and C\ > or 6*6 = C\ = 0). We therefore conclude that 
||ey||r goes to as 5 — > 0. The first three equations of ()135jl implies that e* 5 is a bounded 
sequence in Ho(curl, fi e ) (see definition in the proof of Theorem 16. ip . Therefore, up to 
extracted subsequence, one can assume that e s converges strongly in L 2 (f2 e ) and weakly in 
-£/o(curl, fi e ) to some e € £fo(curl, f2 e ). Passing to the limit as S — > in (|135|) . we observe 
that e G i/o(curl, f2 e ) is solution of 

curl curl e — u> 2 e — 0, in f2 e , 

curl ex n — iLuex xn = 0, on <9f2, 
e x n = on r, 

and therefore e = 0. We then deduce that curl/i* strongly converges to in L 2 (f2 e ). Com- 
ing back to identity (|136p and considering the real part, one deduces after division by 
that 

\\h 5 \\ 2 H(cull ,n e ) < ^ C ||^|U (curI) a e) + C 1 ||curlA tf ||i, (n<) + C 2 ; \Xm (V S ' k h s Tl h s T )\. (142) 

Property l|129jl -(i) shows that 

\lm (V S ' k h s T ,h s T )\ < C 3 5\h s T \l -► 

according to (|14ip . Therefore, considering the limit as S — > in (|142p implies that /i 5 
strongly converges to in £f(curl, O e ), which contradicts H^Hl 2 ^) = 1. □ 
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A Technical Lemmas 

The first Lemma is a slight variation of the classical trace lemma in (curl, O). 

Lemma A.l Let O be a bounded open subset ofM. 3 of class C 2 (and which is locally from 
one side of its normal). Then there exists a constant C depending only on O such that 

||«xn||^_ i(ao) <C|H| i2(0) (||curl«|| La(0) + \\u\\ L , {0) ) Vu e ff(curl.O). (143) 

Remark A.l When O = R 3 + := {y E R 3 / y 3 > 0} , one can easily check, using partial 
Fourier transform in the variables (2/1,2/2); that 

VuG #(curl,R 3 ), ||u x n||^_i < 2 IM^rS) (||curlu|| L2(R 3 } + ||u||i,2( R 3 )). 



) 

Proof. The idea is to see how the proof for the half-space (cf. previous remark) is modified 
when the boundary of O is not flat. 

Let us consider first the cases where there exists a > 0, 6 > 0, S > and h G C 2 (R 2 ) D 
iy 2 <°°(R 2 ) such that 

dOn supp u C E := I ip(y 1 ,y 2 ,0) ; (2/1,2/2) e] - a, a [x] - 6,6 [ j, 
p(]-a,a[x]-&,6[x]0,<y[) C O, 

u is compactly supported in ip^ ] — a, a [x] — 6, 6 [x]0, £ Q , 
where 

^(2/1,2/2,2/3) = (2/1,2/2,2/3 + h(y 1 ,y 2 ))- 
Setting h:=moj3 and n uo (p , one has for 2/3 = 0, 



/(«2 + «3|^-,ii3|^--fii, -fii|^-+fiJU), /:=(1 + |^| 2 )^ (144) 
V aj/2 oyi 02/2 oyi/ 



u x n 



On the other hand, (ei, e2, 63) denoting the canonical basis in the y-sp&ce, 



(Hi) (curlw o ip) ■ e 3 



du 3 


du 3 


dh dii 2 


dy 2 


dy 3 


dy 2 dy 3 ' 


dui 


du 3 


dy 3 dyx ' 


dy 3 


dyi 


du 2 


dui 


9?/3 9/l 

dy 3 dy 2 


dyi 


dy 2 



(145) 
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* - - dh . 
Let us set u 2 = u 2 + u 3 — — , since 



dut du 2 du 3 dh 



d 2 h 



u 3 



9y 3 dy 3 dy 3 dy 2 dy 3 dy 2 
using the formula (|145|) -(i) one gets, setting r := (curlwo ip) = (ri,r 2 ,r 3 ) 



du* 2 
dy 3 



du 3 



d 2 h 



-n 



u 3 



dy 2 dy 3 dy 2 



(146) 



In what follows, we use the Fourier transform T in the variables (2/1,2/2) and denote by 
(6, 6) the dual variable, by Ui the Fourier transform of Ui and u 2 the Fourier transform of 
u 2 . By definition of the norm in i?-5(K 2 ), 



K(v,o)||^_i = / (i + |£IT 5 R(6, 6,0)1^ d6<*6 



Since 



|^(6,6,0)| 2 = -2fte 



ay 3 



(6,6,2/3) ^2/3 



using (|146[) . we have, r\ being the Fourier transform of r\, 



1^(6, 6, 0)| 2 = 2 Tie 



ri «2 



(6,6,2/3) ^2/3 



2 Tie [u 3 u* 2 ] (6 ,6,2/3) */3 

Jo 



2fte 



T(u 3 



d 2 h 
dy 3 dyi 



(6,6,2/3) ^2/3 



We divide the above equality by (1 + 1£ | 2 ) 2 and integrate over £. Next we use (1 + |£| 2 ) 2 < 
1, 161 (1 + l?| 2 )~^ < 1 and Plancherel's theorem to obtain, since h € M^ 2 '°°(M 2 ), 



||^( V) 0)||^_i <2 



Js. 2 



\n\ \u* 2 \dy + C 



\u 3 \ \u* 2 \ dy. 



Coming back to the variable x through the change of variable x = f(y), we easily get, since 
K| < \u 2 \+C\u 3 \ 

IKO.-. 0)11^-4 < c (\Hh(o) + \\u\\l2(0)\\cuiIu\\ L 2 (0 )) 
where the constant C only depends on h. 
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Finally, using the lemma YK72\ (notice that / := (l + |-D^| 2 ) 2 belongs to W 1,oc ), we get 



11/(^ + ^3^)11^^ < C (\\u\\ 2 L . 2{0) + |MU 2(0 )||curhi|U 2(0) ) (147) 
In the same way, one obtains 

||/(fi 3 ^.-u 1 )|| H _ i < c(|| U ||| 2(0) + || U || L2(0) ||curl U || i2(0) ) 

B h (W8) 

H/(-"i^+«2^)|| H -i (r) < c(||u||2 2(0) + || u || L2(0) ||curlu|| L 2 (0) ) 

Thanks to (|144|) and by definition of the norm in H~*(T), estimates (|147j) and (|148|) lead 
to the desired inequality. 

Obtaining the same inequality in the general case can be deduced by the using a partition 
of unity (<Pi)i=i n of O and noticing that 

||curl<piu||.£3(o) = \\<Pi curlu + Vtpi x u||ia( ) < \\fi\\oo ||curMU»(o) + ||Vy>i||oo IMU 2 (o)- 

□ 

Lemma A.2 Let f E W 1 '°°(R n ) and g E H-^(R n ) then f g € H-^(R n ) and one has, 

\\fg\\ H - i(Rn) <^ \\g\\ H . i(Rn) . 

Proof. We first remark that if -0 € Hi(R n ), then 

\\fip\\ Hi{Rn) <^ wfiiw^) \m Hi{Rn) 

which is deduced by interpolation from the (obvious) inequalities: 

' V0 e L 2 (R"), \\f^\\ H i (Rn) < IMUw, 

Next, if g € H-^(R n ), we have, for any ip E H?(R n ), 

\(fg,i>)\ = \(g,M)\<\\g\\ H - i(Rn) \\fi>\\ ffiiRn) <3i \\g\\ H - i(Rn) 11/11^-^) M H - iiRn) , 

from which one easily concludes. □ 
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Lemma A. 3 Let O C R 3 be a bounded open set with a C 2 boundary T. There exists a 
constant C that depends only on V such that 

\H H h {r) < C ^r u W H -h { r) + H U W)) v w e ff*(r). 

Proof. In the case O — {(xi, X2, £3) £ R 3 / X3 > 0} one can check by using Fourier transform 
in the plane {x\,x<x) that 

||u|| 2 ! = ||V r w|| 2 1 + IMI^m ■ 

The inequality is therefore trivially verified in this case. The general case can be easily 
deduced by using local parameterizations of the boundary T. This is where the C 2 -regularity 
of r is taken into account. □ 

Our next lemma is a sharper version of classical compact-embedding theorem for spaces of 
L 2 functions with bounded divergence and curl into L 2 . We set 

ff(curl, div, O) := {u G L 2 {Of / curlu G L 2 {Of and div u G L 2 (0)} 

equipped with the norm 

IMlH(curi,div,o) = \\ u \\h(0) + H c urH| 2 2(0) + ||divu|| 2 2(0) . 

Lemma A. 4 Let O C i? 3 be a bounded simply connected open set with C 2 boundary T. 
Then every bounded sequence (uk)keft of H (curl, div, O) such that 

_ 1 

(ufe|r x n)ken is convergent in H t 2 (r) (149) 
has a convergent subsequence (uk>) in L 2 (0) 3 . 

Proof. Our proof is an adaptation of the proof given by Costabel in the case where, instead 
of (|149[) . one has an L 2 control of the boundary term of the sequence (see Theorem 2 of [S]). 

The idea is to make a Helmholtz decomposition of Uk of the form: 

Uk=w k + Vpk, (wk, Pk ) G H^O) x L 2 (0), divw fe = 0, (150) 
constructed in such a way that: 

(i) Wk is bounded in H 1 (0) (and thus admits a converging subsequence in L 2 (0) 3 ) : this 
uses the fact that curl-it^ is bounded in L 2 (0) 3 , 

(i) Vpfc admits a converging subsequence in L 2 (0) 3 : this uses the fact that divitfc is 
bounded in L 2 (0) 3 and that (uk x n)| r converges according to (|149[) . 
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In order to construct Wk from curl Uk we first construct an extension of curl Uk in R 3 which 
has a compact support (independent of k) and is divergence free. For this, we choose a ball 
B containing O in its interior and will constitute the support of the extension of curl Uk- 

First notice that since curl it^ G ff(div; O), the trace curlu^ • n\ r is well defined in H~^(T) 
and satisfies, since divcurlufe = 0: 

||curlu fc • n\\ H _i < C ||curlu fe || £ 2( 0) , (curlu fe • n, l) r = 0. (151) 

Therefore, the following Neumann problem in B\0: 

Aip k = in B \ O, 

dnfk = (curlwfe) • n on 30, (152) 
k d n if k = on dB, 

admits a solution ipk £ H 1 (B \ O), which is unique if we impose in addition 

tpk dx = 0. 



'B\0 

Moreover, using the Poincare-Wirtinger inequality and (|151|) 

IV^/cl 2 dx = ( (curlufc) • n,lpk~) r < C ||curlw fc || L 2 (0 ) \\Vtpk\\ L 2( B \o) 

b\o 

from which we deduce 

ll v ^|| L 2 (B \o) ^ C \\ CUTlu k\\Li(o)- (153) 
We then introduce the extension \k of curl Uk as 

in O, 

Xk = { Vtfk in B\0, 

in R 3 \ B. 

Of course, \k € L 2 (R 3 ), is compactly supported in B and satisfies, thanks to p53|) 

||Xfc|U 2 (R 3 ) < C ||curlw fc || L 2 (0) . 

Moreover, the two boundary conditions in (|152p have been chosen in order to enforce the 
continuity of Xk ■ n across dO and dB so that: 

divxfc^O inR 3 . 
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Next, we introduce as the unique solution of he following Laplace problem in IR 3 

-AV k =Xk inM 3 , * fe € Hf oc (R 3 ) 3 , Vf,el 2 (l 3 ) 3 . (154) 
It is well known that is given by: 

= G * Xfe, G = — -r— r, (the fundamental solution of the Laplace operator) 
47r|a;| 

and satisfies in particular 

ll*fe||ffi(o) < c WxkW^m • (155) 

Moreover, 

div X fe = div*fe = 0, (156) 

Next we define 

w k =cw\y k , e L 2 (R 3 ) (157) 
whose restriction to O is the good candidate for (|150|) . Indeed 

curl Wfe = curl (curl* fc ) = grad (div* fe ) - A* fe = Xk in M 3 , (see (|154|) and (|156|l ) 

which implies in particular 

curl Wfe = curl u k in O. 

Before constructing p k , we first check property (i). The fact that Wk is bounded in L 2 (0) 3 
results directly from (|155[) . Next, we show that w k is bounded in i? 1 (R 3 ), which implies in 
particular (i). Indeed using the Fourier transform in M. 3 we deduce from (|157[) and (|154[) 
that (£ denotes the dual variable of x and u the Fourier transform of u) : 

\e imoi 2 = '^i ' 2 < i»(oi 2 

which yields, by Plancherel's theorem 

llVlWfcllla^S) < ||Xfc|||3(R8) < C ||curlUfc||| 2(0) . 

From now on, we can therefore assume that (up to extracted subsequence) w k converges in 

Since curl(ttfc — w k ) — and O is simply connected, one can construct p k (unique up to an 
additive constant) such that Vpfc = u k — w k (use for instance Theorem 2.9 of [1]). Fixing 
Pk by imposing that J p k dx = gives raise to a bounded sequence p k in 7J 1 (0) by the 
Poincare-Wirtinger inequality. Since we further have that (div it/.) is bounded in L 2 (0) and 
(u>k) is bounded in H (O), then, up to extracted subsequence, one can assume that divufe is 
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convergent in H 1 (0), Wk\do is convergent in H 2 (dO) and pk\ao is convergent in L 2 (dO). 
We shall deduce that pk is strongly convergent in We first observe that pk satisfies 

-Ap k = divu fe , in O, 

V^fc x n — Uk x n — Wk X n, on 90. 
Let m and A; be two indexes. From 

A(p fc - p m ) = div(u fe - u m ) in O, 

(158) 

V(pfc - p m ) x n = (u k - u m ) x n - (iu fe - w m ) x n on dO. 

and using the classical theory for elliptic equations one gets the existence of a constant C\ 
such that 

II Vp fc - Vp m || L 2 (n) < Ci M|div(it fc -u m )\\ H -i( ) + \\Pk ~Pm\\ H ^, d0) )■ (159) 
On the other hand, using Lemma lA.3l one has 

\\Pk -Pm|| H i (ao) < C2 (||V(p fc -p m ) x njl^i + |bfe -p ro |U=(90)J • (160) 
From the second equation of (|158[) , (|159[) and (|160p it is easily seen that 

IjVpfc - Vp m || L 2 (0) < C 3 (j|div(u fe - u m )|| ff -i ( o) + \\(uk - u m ) x n\\ H _i 

+ \\w k -W m \\ H -^ igo) + \\Pk -PrnWmdO)) ■ 

Using assumption (|149|) one concludes Vpk is a Cauchy sequence in L 2 (dO). The result of 
the lemma is then proved since Uk — Wk + Vpfc. □ 

Lemma lA. 41 also applies to domains Qi that are not simply connected. This is proved in the 
following lemma. 

Lemma A. 5 The result of Lemma \A.4\ applies to bounded open domains O € R 3 of class 
C 2 . 

Proof. Let x be an arbitrary point in O. If x G O, one defines U x as a ball centered at x 
such that U x C O. If not, one defines U x as a neighborhood of x such that there exits a 
bijective map 4> x Q <-^> U x such that 

^eC 1 ©, f^c 1 ^), W+) = £/xno, et (i>(Q ) = undO, 

where Q denotes the unit cube of R 3 , Q+ := {x 6 Q | x 3 > 0}, and Qo = S Q | X3 = 0}. 

With this definition one observes that U x DO is a simply connected domain for all x € O. 
By the compactness of O one can extract a finite covering of O from {U x ;x G O}. Let us 
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denote by {[/,-, i£ 1} this finite covering and consider a partition of unity (9i) ie j C C°°(R 3 ) 
subordinated to this covering, i.e. 

supp 6i C Ui , $i = 1 on O. 

iei 

Then define u l n := 9i u n for all i 6 I. It is easy to see that for every i, the sequence 
(u l n ) satisfies the hypotheses of Lemma fA.41 with O replaced by Ui. Using a finite diagonal 
process, one can therefore assume that there exists a subsequence such that 

u\ k converges in L 2 (Ui) for all i e I. 

Consequently, the sequence u Uk = J2iei u li k ^ s convergent in L 2 {0). □ 
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